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A  method  for  predicting  swell  from  tropical  cyclones 
using  a  spectral  wave  model  (TYWAVES)  was  tested.   The  model 
was  applied  to  predicting  swell  propagating  from  three 
typhoons  in  the  Western  North  Pacific  through  gaps  in  the 
Ryukyu  Islands  into  a  region  of  the  East  China  Sea.   The 
model  involves  a  source  region  concept  which  considers 
only  the  swell  emanating  from  regions  of  peak  energy  in 
moving  typhoons.   For  three  representative  typhoons, 
predicted  heights  were  not  significantly  different  from 
the  observed  heights.   The  time  of  occurrence  of  the  pre- 
dicted peak  height  agreed  well  with  observational  values 
for  the  swell  from  two  typhoons,  but  lagged  by  6-12  hours 
for  the  third. 

The  dominant  swell  period  and  direction  predicted 
by  the  model  were  not  verifiable  by  data  available  for 
this  study. 

Shoaling  and  refraction  effects  were  considered  in  the 
prediction,  in  a  simplified  way,  but  attenuation  was 
ignored  even  for  the  passage  of  energy  through  the  Ryukyu 
Islands. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 11 

A.  OBJECTIVE  OF  THE  STUDY 11 

B.  TYPHOONS  IN  THE  WESTERN  NORTH  PACIFIC 11 

C.  TYWAVES  MODEL  (THE  TYPHOON  WAVES  PROGRAM) 12 

D.  FORECASTING  AND  VERIFICATION 13 

E.  TERMINOLOGY 22 

II.  PREDICTION  OF  TROPICAL  STORM  WAVES  BY  TYWAVES 
MODEL 2  5 

A.  LOCATION  OF  THE  SELECTED  POINT-SOURCES 25 

B.  SEA-SPECTRA  AT  THE  POINT-SOURCES 2  8 

C.  SELECTION  OF  THE  DIRECTIONAL  ENERGY  AT 
POINT-SOURCES 2  9 

D.  TOTAL  AND  DIRECTIONAL  ENERGY  AT  POINT- 
SOURCES 3  0 

1.  Typhoon  Hope 31 

2.  Typhoon  Irving 31 

3.  Typhoon  Owen 32 

III.  PREDICTION  OF  SWELL  CHARACTERISTICS 36 

A.  THE  SETTING  OF  THE  FORECAST  SITE  AND  ITS 
RELATION  TO  THE  WAVE  SOURCES 36 

B.  GROUP  VELOCITY 37 

C.  TRAVEL  DISTANCE  AND  ARRIVAL  TIME 3  9 

IV.  THE  ENERGY  PROPAGATION  PROCEDURE 45 

A.   PROPAGATION  OF  SPECTRAL  ENERGY  COMPONENTS 45 


B.  SHOALING  AND  REFRACTION  FACTOR 4  7 

1.  Shoaling  Factor 47 

2.  Refraction  Factor 50 

C.  SHOALING  AND  REFRACTION  OF  THE  SPECTRAL 
COMPONENTS 51 

1.  Shoaling  Process  (Computations) 52 

2.  Shoaling  Energy  Components  From  Each 
Source  Versus  Arrival  Time 52 

3.  The  Predicted  Swell  Waves  at  Prediction 
Site 5  3 

D.  OBSERVED  DATA 53 

E.  COMPARISON  OF  SWELL  PREDICTIONS  AND 
OBSERVATIONS 6  3 

F.  ERROR  SOURCES 64 

V.     CONCLUSIONS 67 

APPENDIX  A  (Outputs  of  TYWAVES  for  Typhoon  Owen  on 

12  GMT  September  26) 68 

APPENDIX  B  (Sea  Spectra  at  Source  Region  for  Three 

Typhoons  and  Significant  Wave  Distribution)-  76 

APPENDIX  C  (Propagation  Energy  Arrival  Time) 89 

APPENDIX  D  (Shoaling  of  the  Spectral  Energy  Components)  -  94 

LIST  OF  REFERENCES 113 

INITIAL  DISTRIBUTION  LIST 115 


LIST  OF  TABLES 

1.  1979  Significant  Tropical  Cyclones 14 

2.  1979  Significant  Tropical  Cyclone  Statistics 14 

3.  Best  Track  for  Super  Typhoon  Hope 19 

4.  Best  Track  for  Typhoon  Irving 20 

5.  Best  Track  for  Typhoon  Owen 20 

6.  Typhoon  Centers  and  Selected  Points  Location 2  7 

7.  Period-Directional  Spectrum  at  Point  4  of  06  GMT 
September  25,  1979 30 

8.  The  Visual  Observation  Data  by  NODC  and  ROK  Navy- 6  0 


LIST  OF  FIGURES 

1.  Typhoon   Tracks 15 

2.  Typhoon  Hope,  Best  Track 16 

3.  Typhoon  Irving,  Best  Track 18 

4.  Typhoon  Owen,  Best  Track 17 

5.  Possible  Refraction  Diagram  at  Point  4  at  06  GMT 
September  25,  1979 29 

6.  Typhoon  Hope,  Total  and  Directional  Energies  at 
Point  Sources 33 

7.  Typhoon  Irving,  Total  and  Directional  Energies  at 
Point  Sources 34 

8.  Typhoon  Owen,  Total  and  Directional  Energies  at 
Point  Sources 35 

9.  Typhoon  Hope  Swell  Arrival  Time 42 

10.  Typhoon  Irving  Swell  Arrival  Time 43 

11.  Typhoon  Owen  Swell  Arrival  Time ■ 44 

12.  Typhoon  Hope,  Directional  Energy  Propagation 54 

13.  Typhoon  Hope,  Total  Propagated  Energy  Comparison 5  7 

14.  Typhoon  Irving,  Directional  Energy  Propagation 55 

15.  Typhoon  Irving,  Total  Propagated  Energy  Comparison-  58 

16.  Typhoon  Owen,  Directional  Energy  Propagation 56 

17.  Typhoon  Owen,  Total  Propagated  Energy  Comparison 59 


LIST  OF  MAP  AND  PLATE 


MAP  1.   East  China  Sea  Bottom  Contours 21 

PLATE  C-l.   Illustration  of  Various  Functions  of  =—  4  8 

o 

FIGURE  2-19.   Changes  in  Wave  Direction  and  Height  Due 
to  Refraction  on  Slopes  with  Straight 
Parallel  Depth  Contours 49 


ACKNOWLEDGMENT  S 

I  wish  to  express  my  deepest  thanks  to  Professor  J.  B. 
Wickham  for  the  concepts  and  methods  in  this  thesis.   The 
completion  of  this  work  is  due  to  his  effort,  experience 
and  guidance. 

I  also  wish  to  recognize  the  invaluable  cooperation 
and  assistance  of  Samson  Brand,  Naval  Environmental 
Prediction  Research  Facility,  Monterey,  California,  9  3940; 
Kevin  M.  Rabe,  Science  Applications  International,  Monterey, 
California,  9  3940;  and  of  Professor  J.  Von  Schwind,  Naval 
Postgraduate  School,  Monterey,  California,  9  3940. 


10 


I.   INTRODUCTION 

A.  OBJECTIVE  OF  THE  STUDY 

A  serious  problem  for  naval  activities,  ocean  industry, 
shore  protection  and  fisheries  in  the  Korean  south  coastal 
area  is  the  lack  of  adequate  estimation  of  the  waves  re- 
sulting from  typhoons.   The  objective  of  this  study  is  to 
make  a  partial  test  of  the  TYWAVES  for  forecasting  swell 
from  tropical  cyclones  which  arrive  at  a  single  observa- 
tion site  near  Cheju-do,  Korea  (33.2°N  126. 6°E). 

Forecasting  based  on  wave  fields  predicted  in  a 
typhoon  area  by  the  TYWAVES  [13]  developed  by  NEPRF 
(Naval  Environmental  Prediction  Research  Facility)  are  to 
be  verified  against  swell  observations  made  near  Cheju-do. 

B.  TYPHOONS  IN  THE  WESTERN  NORTH  PACIFIC 

During  1979,  the  Western  North  Pacific  experienced 
2  8  tropical  cyclones.   Table  I,  from  "1979  Annual  Typhoon 
Report"  [2],  shows  the  significant  tropical  cyclones  for 
that  year.   Table  II  from  [2]  shows  the  monthly  distribu- 
tion of  tropical  cyclones  in  1979  and  other  statistics. 

Most  typhoons  occur  in  the  summer  season  and  their 
tracks  can  be  classified  as  one  of  three  typical  tracks. 
The  first  type  is  that  passing  south  of  Taiwan  toward  the 
west,  the  second  is  that  crossing  over  Korea  through  the 
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Ryukyu  Islands,  and  the  third  is  that  passing  east  of 
Ryukyu  Islands. 

I  selected  one  typhoon  of  each  type;  one  each  in  July, 
August,  and  September  in  1979.   They  are  Typhoons  Hope, 
Irving  and  Owen  shown  in  Figures  1,  2,  3  and  4  from  [2]. 
tropical  cyclone  "best  track"  information  is  shown  in 
Tables  III,  IV,  and  V  for  Typhoon  Hope,  Irving  and  Owen, 
respectively. 

C.   TYWAVES  MODEL  (THE  TYPHOON  WAVES  PROGRAM) 

The  SOWM  (Spectral  Ocean  Wave  Model)  run  at  FNOC 
(Fleet  Numerical  Oceanography  Center)  utilizes  a  coarse 
operating  grid  system  (100~190NM)  which  does  not  allow 
sufficient  resolution  to  describe  adequately  the  resultant 
sea  state  in  typhoon  areas.   Thus,  TYWAVES,  an  improved 
model  for  typhoons  with  a  locally  finer  grid,  was  developed 
by  NEPRF.   FNOC  judged  that  TYWAVES  is  worthy  of  evalua- 
tion as  a  possible  operational  typhoon  sea  state  model 
[14]. 

The  TYWAVES  is  intended  to  produce  fields  of  signifi- 
cant wave  height  and  spectral  wave  properties  on  a  mesh 
size  consistent  with  the  scales  of  tropical  cyclones  and 
is  designed,  primarily,  for  the  application  to  the 
western  North  Pacific. 

The  detailed  outputs  of  TYWAVES  are  in  the  form  of 
fields,  at  each  of  12  points,  of  spectral  energy 
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components,  significant  wave  heights,  maximum  wave  periods, 
and  the  predominant  wave  directions,  for  00,  12,  24,  48  and 
72  hours,  where  00  hours  is  the  time  of  the  first  typhoon 
warning  issued.   An  example  of  those  outputs  for  Typhoon 
Owen  is  shown  in  Appendix  A. 

Using  these  spectral  energy  components  at  selected 
source  points,  I  was  able  to  make  predictions  based  on  the 
propagation  of  these  components  as  swell  into  the  region, 
south  of  Cheju-do,  Korea  (see  Map  I,  page  21). 
More  details  of  TYWAVES  are  described  in  Ref.  11. 

D.   FORECASTING  AND  VERIFICATION 

Extremely  high  sea  states  are  known  to  be  generated  in 
the  quadrant  to  the  right  of  the  direction  of  movement  of 
typhoons.   The  wave  generation  in  that  quadrant  of  the 
typhoon  derived  in  TYWAVES  from  a  spectral  model  utilizing 
the  Pierson-Moskowitz  (1964)  spectrum.   The  model  describes 
the  spectroangular  components  of  the  waves  present  at  a 
number  of  grid  points  in  the  region  of  strong  winds.   Each 
spectral  component  of  interest  is  then  permitted  to 
propagate  at  its  appropriate  group  velocity  to  the  fore- 
cast site.   The  method  is  applied  to  three  western  North 
Pacific  typhoons  in  1979  and  the  forecast  products  are 
compared  with  the  observed  swell  data  from  the  south 
coastal  area  of  Cheju-do,  Korea. 
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TAB>-t   1 .                             WESTERN 

NORTH  PACIFIC 

19/9  SIGNIFICANT  TROPICAL  CYCLONES 

! 

CALENDAR 

MAX 

MIN 

NUMBER 

PERIOD 

DAYS  OF 

SFC 

08S 

OF 

DISTANCE 

CYCLONE     TYPE     NAME        OF  WARNING 

WARNING 

WIND 

SLP 

WARNINGS 

TRAVELLED 

01         TY      ALICE      01  JAN-14  JAN 

14 

110 

930 

51 

2597 

02         TY      BESS       20  MAR-25  MAR 

6 

90 

958 

21 

1804 

03        TY      CECIL      11  APR-20  APR 

10 

80 

965 

40 

2535 

04        TS      DOT       10  MAY-16  MAY 

7 

40 

984 

24 

2876 

05        TD      TD-05      23  MAY-24  MAY 

2 

30 

998 

6 

2170 

06        TY      ELLIS      01  JUL-06  JUL 

6 

85 

955 

22 

1612 

07        TS      FAYE      01  JUL-06  JUL 

6 

40 

998 

20 

1837 

08        TD      TD-08      24  JUL-25  JUL 

2 

20 

1004 

5 

1264 

09         ST      HOPE       27  JUL-03  AUG 

10 

130 

898 

33 

3928 

10         TS      GORDON     26  JUL-29  JUL 

4 

60 

980 

13 

1058 

11         TO      TD- 1 1      03  AUG-06  AUG 

4 

25 

997 

14 

1088 

12         TY      IRVING     09  AUG-18  AUG 

10 

90 

954 

38 

2732 

13        ST      JUDY       16  AUG-26  AUG 

11 

135 

887 

39 

2502 

14         TD      TD-14      18  AUG-20  AUG 

3 

20 

1006 

9 

605 

15        TS      KEN       01  SEP-04  SEP 

5 

60 

985 

13 

1418 

16         TY      LOLA       02  SEP-08  SEP 

7 

90 

950 

23 

1298 

17        TY      MAC       15  SEP-24  SEP 

10 

70 

984 

35 

1831 

18         TS      NANCY      19  SEP-22  SEP 

4 

45 

993 

14 

528 

19        TY      OWEN      22  SEP-01  OCT 

10 

110 

918 

37 

2151 

20        TS      PAMELA     25  SEP-26  SEP 

3 

45 

1002 

6 

984 

21        TS      ROGER      03  OCT-07  OCT 

6 

45 

985 

16 

1920 

22        TY      SARAH      04  OCT- 15  OCT 

12 

no 

929 

43 

1194 

23        ST      TIP       05  OCT-19  OCT 

16 

165 

870 

60 

3972 

24        ST      VERA      02  N0V-07  NOV 

6 

140 

915 

23 

1868 

25        TS      WAfNE      08  NOV- 13  NOV 

6 

50 

990 

22 

1559 

26        TD      TU  26      01  OEC-02  DEC 

2 

30 

998 

6 

1070 

27         TY      ABBY       01  DEC-14  DEC 

14 

110 

961 

52 

4044 

28        TS     bEN       21  DEC-23  DEC 

J 

60 

990 

10 

2245 

1979  TOTALS 

149* 

695 

•OVERLAPPING  DAYS,  INCLUDED  ONLY  ONCE  IN  SUM. 

TABLE   2. 

J 

1979  SIGNIFICANT  TROPICAL  CYCLONE 

STATISTICS 

| 

'8)1 

WESTERN 

(1959- 

NORTH  PACIFIC      JAN   FEB   MAR   APR   MAY   JUN 

JUL   AUG 

SEP 

OCT   NOV 

DEC 

TOTAL     AVERAGE 

I 

TROPICAL 

j 

DEPRESSIONS        0O0O1O 

1     2 

0 

0     0 

1 

5         4.8 

* 

TROPICAL  STORMS     0     0     0     0     10 

2     0 

4 

1     1 

1 

10        10.0 

i 

TYPHOONS           10     1     10     0 

?           2 

? 

2     1 

1 

13       18.0 

i 

ALL  CYCLONES       10     112     0 

5     4 

6 

3     2 

3 

28       32.8 

• 

(1 

959-78)  AVERAGE   0.6   0.4   0.6   0.9   1.4   2.1 

5.2   6.8 

6.0 

4.8   2.7 

1.3 

32.8 

FORMATION  ALERTS    23  of  the  27  (85%)  Formation  Alert  Events  developed  Into  tropical  cyclones. 
5  of  the  28  (18t)  tropical  cyclones  did  not  have  a  Formation  Alert. 


WARNINGS 


Number  of  warning  (lays:   149 

Number  of  warning  days  with  ?  cyclones:   38 

Number  of  warning  days  with  3  or  more  cyclones:   S 
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Figure!  Typhoon    tracks 
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Figure  3- 
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IT      f'RST   WARNING   'SSUEO 
-.1      LAiT    WARNING     SSUEO 
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Table  3.    Best  track 


SUPER  TYPHOON  HOPE 


90/J4/MJ 

0'2*l2t 

oUnUJ 

0 7231B/ 

o  '23002 
072BU6/ 
u7<:3i2/ 

aii!  i  mi 

0  l£l  <jbl 
U72/122 
j'Ulnl 

0  729U42 

J729I22 
0 '2911/ 

j  72.062 
072*122 
Ol<Hl*i 

j  ' i  1001 

,'jju'; 

07io 182 
4'Jlo02 

li'JUii 

ill  H  ttil 

)8U  U3i 
4»0 1 1^/ 

48o  i I  a/ 

J842122 

464)00/ 
4B4i06/ 
090il22 

oau Jia/ 
<jau»uu/ 
o8o»o6.; 
084.12/ 
ubo. ia/ 

o803j6/ 

Jlu)l.'.' 
0BO31B/ 

4846002 

0B0646/ 

.'■on; 

UBO 704/ 

osorizj 

jbo  I  ,  a/ 
oaoauO/ 
jaoaob/ 
uao9i2z 


S«si    Tj|^ 

H  .IV. 

1*7. a  /I 

146.9  <n 

I  «  6  .  ?  i  A 

14S.S  2. 1 

I44.ii  iS 


l  4  I 


i 


1,2.4 
141. A 
I  40.  I 
1  19.9 
1  <.0  .  I 
140.1 

im.k 

t   18. A 

I  17.  a 
1  17. s 
I  16.9 
I  16.3 
I  IS.  7 

I    IS. 4 

1  14. a 

114.? 
1  13. S 
1  12.7 
111.4 
1  10.* 
I  '"».» 
I  >7.a 
l->6.  / 

l?4.7 

l?3.? 

Ml.* 
I'd.  I 
I  14.? 

116.4 

lll.i 
111. 7 

119.* 
1"7.S 
1  IS. 4 
II).  I 
1 .1 1  .  7 

111.4 

19.? 

•  7.9 
•>6.7 

<S.» 
94.9 
93. S 

>2.7 
92. A 

>1.4 

'O.-i 
90.1 
9*.  7 
»9.n 
18.1 
97.? 

96.4 

iS.s 
4*.* 


I  IS 
I  3i 

l  tn 
l  in 
|2s 
121 
US 

ll'S 

■>s 

71 
*,l 
IS 
11 
?1 


?1S 

U.I 

0.0 
O.n 
o.o 
o.n 

11. n 

li.i 

lt.i 

11.4 
II.  1 

12.0 

12.7 

I  i.  7 

0.0 

o.o 
la. 2 

1  7.5 

14.2 
19.1 
16.9 
16.2 
16.1 
16.9 
la.  4 
'.7,1 
17.2 
17.5 
19. S 
14.2 
l».7 
20.1 
20.7 

20.  M 
21.4 

21.9 
22. n 
22.4 

22. b 
22.7 
22.  5 
0.0 
0.1 
0.0 
9.0 
0.1) 
0.3 
0.0 
0.0 
0.0 
4. ft 
o.n 

0.0 

o.o 
o.n 

0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.1 
0.0 


.««    ,1- 

.INI 


24     M<1i 


0. 

11 

u. 

0 

1 

0. 

0. 

a 

u. 

0. 

n 

0. 

.« 

1 

'b. 

41. 

1 

2*. 

»2 

7 

2D. 

41  . 

s 

2u. 

40 

A 

20. 

)■) 

7 

2u. 

19 

7 

2". 

40 

7 

20. 

0< 

9 

0. 

Oi 

1 

4, 

17 

•a 

2b. 

3H. 

9 

2$. 

37 

? 

2b. 

]« 

. 

2b. 

lb 

? 

Jb. 

3S. 

1 

*". 

3. 

9 

65. 

3. 

4 

'u. 

31 

* 

'b. 

3? 

* 

10. 

32 

1 

»u. 

31 

1 

*b. 

!* 

1 

I  40. 

iH 

A 

n'b. 

?6 

9 

I   3j. 

>.% 

* 

I  JO. 

J3 

? 

I  ill. 

^l 

S 

'2b. 

20 

1 

i2u. 

ia 

4 

1  13. 

16 

s 

i  ob. 

14 

1 

•»u. 

12 

1 

'o. 

10 

I 

60. 

OH 

A 

ib. 

0 

1 

0, 

0 

A 

0. 

0 

A 

u. 

0 

A 

u. 

1 

A 

o. 

0 

A 

u. 

0 

A 

0. 

0 

1 

0. 

0 

1 

0. 

0 

1 

1). 

0 

A 

J. 

0 

1 

u. 

0 

1 

Ii. 

0 

1 

o. 

0 

>A 

0. 

0 

►  1 

i. 

0 

1 

0  . 

0 

■  n 

0. 

ft 

|A 

4. 

0 

■  A 

ft. 

0 

ii 

0. 

0 

1 

0. 

1  3. 
11. 


N14& 
•  IN) 


o  . 
0. 


1«. 
1  i. 


19. 
-ft. 
-0. 

S. 
B5. 
62. 
Ill, 
29. 
29. 

a. 
.«, 

a. 

9. 

1  7. 


0. 

0. 


-10. 

-it). 


5. 

il. 
-lo. 


9. 

la. 
1 1. 

9. 
21. 

29. 
13. 

-ft. 
-I). 
-0. 
-0. 
-0. 


-ft. 

-0. 

-ft. 


-0. 
-0. 

-0. 


ID. 


0. 
0. 


PnSir 

0.1  fl 

o.n       n 

0.1 

ft.n 


H     fOHtt.* 


O.i 
O.i 
0.1 

lH.i 

?n.  7 

21.* 
?2.? 

I  7.S 
16.? 
17.1 


I  31. 

13*. 

13». 

131. 

13?. 

13?. 

13?. 
17.1  131. 
18.1    I3i. 


|8.3 
18.1 

19.1 

>O.S 

21.7 

'l  .* 

?2.A 

?2.* 
?2.7 

>  3.1 

?3.4 

?2.» 
?3.i 
O.n 
O.n 
O.n 
O.i 
O.i 
O.n 
0." 
O.n 
O.n 
O.n 
O.n 
O.n 
O.n 
O.n 
O.n 
O.n 
O.n 
O.n 
O.n 
O.n 
O.i 
O.n 
O.n 
O.n 
O.n 


12«. 

129. 

12". 

124. 

121. 
I2n. 
111. 
117. 
II*. 
114. 
I  1>. 
I  H. 
I  0*. 


■  I  911 
0. 
0. 
0. 
0. 


0 

n 

n  .  II 

0. 

12 

s 

I4n.» 

JO. 

|2 

i 

1  3".« 

Jft, 

|2 

? 

1  3*.  1 

JO. 

I  2 

7 

1  37. M 

JO. 

|3 

n 

1  3  7.11 

)0. 

13 

* 

13*.  7 

JO. 

I  4 

. 

|3S.» 

<;i. 

Oil 


)o 

Id 

171 

14< 

I.. 

1  u 


h»Kii 


'2    iJUH    fn9»C4«l 


0. 

0. 

0. 
•0. 
•0.    i 
ii.    J 
J5.    JJ 


bO. 

an. 

75. 

»0. 

90, 

90. 
I  10. 
I  10. 
110. 
I  10. 
120. 
1/0. 

90. 

'a. 


-u 


1 1>0 


II) 


121 

IV, 

1  ,/ 


29 
5 II 

116 

m.    lbs 
25.    223 

•a.   ibj 


JO 
0. 
0. 

o. 
o. 
o. 

0. 
0. 
0. 
9. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 


Ill) 


9i>si  r 

n.il  0 

n.n  0 

n.il  0 

n.il  0 

n.n  0 

II."  I J7 

II.*  I J4 

11.*  I J» 

11.9  1 JJ 

14.1  I J2 

14.9  | j) 

I*,  il  1J1 

n.il  0 

ii .  il  0 

n .  ii  0 

I  9.  '  129 
133 


>1.  I 


ii    131 


J».' 


9.9 


129 

4  124 
I  1  Jt) 
I    129 

I  128 
7     126 

126 
123 

?n.?    124 

-".'i    120 

•".9  117 

II  113 
il  lib 

,.  0 

'I  0 


?4 


•  0 
.0 

0. 
0. 

.0 

ft. 

•  1) 

4. 

•  0 

J. 

•  0 

J3. 

.  * 

J3. 

.d 

Ja. 

•  » 

i3. 

.  7 

i3. 

.  4 

33. 

.  4 

30. 

.4 

4. 

.0 

1). 

•  u 

4. 

.  . 

3U, 

•  1 

34. 

.  J 

43. 

.a 

•3. 

•  0 

60. 

•  1 

64. 

•  J 

13. 

.  1 

1). 

•  i 

1  40. 

.  i 

100. 

.8 

|2o. 

■  . 

120. 

.0 

1  40. 

.4 

ao. 

•  / 

3a. 

•  1 

23. 

.4 

4. 

.J 

J. 

.0 

0. 

.4 

0. 

.J 

4. 

■  4 

I). 

.4 

4. 

•  0 

). 

•  0 

J. 

•  0 

0. 

.0 

I)  . 

•  9 

4. 

.0 

0. 

.4 

u  . 

.4 

4. 

.0 

4. 

•  4 

4. 

•  4 

4. 

•  0 

0. 

•  4 

4. 

.0 

0. 

•  0 

0. 

.j 

0. 

.0 

4. 

•  0 

0. 

•  1) 

0. 

.0 

4. 

.1) 

4. 

.0 

0. 

.0 

4. 

.4 

0. 

Ijl 

•  1  NO 

was 

1 1 

-4. 

0. 

1.0 

0.0 

-4. 

0. 

1.0 

O.'l 

-0. 

0. 

n.O 

0.4 

-0. 

ft. 

n.  a 

0.4 

-4. 

ft. 

1.4 

4.0 

1  *7. 

15. 

IS.. 

I  H.2 

2*6. 

10. 

1*..4 

I  I0.3 

•"■*. 

5. 

14.8 

I  14.4 

>9*. 

0. 

19.2 

i'9.<; 

3114. 

0. 

19.7 

1^9.3 

313. 

0. 

1*.7 

l?7.U 

2S4. 

•10. 

17.8 

127.3 

-0. 

ft. 

n.O 

0.0 

-0, 

0. 

n.O 

4.0 

-4. 

0. 

n.O 

0.4 

290. 

-2b. 

21.3 

1/3.2 

!*■». 

•30. 

24.} 

123.6 

193. 

•40. 

2*. 8 

l?7.2 

197. 

•*b. 

2*. 8 

1/3.3 

46. 

•  4ft. 

19.4 

124.  1 

I«3. 

-as. 

19.6 

127.2 

1«2. 

•45. 

2n.3 

1/3.7 

240. 

•*S. 

21.6 

1/3.1 

21  4. 

•JO. 

2n.» 

1/2.1 

?*». 

•25. 

21.0 

1/2.3 

il4. 

0. 

21  .0 

1/1.8 

194. 

5. 

?i.3 

1/4.9 

214. 

-5. 

2L0 

lla.2 

2'2. 

-5. 

21.8 

lli.8 

211. 

•JS. 

1.0 

0.0 

113. 

•25. 

i.O 

0.0 

-0. 

0. 

n.O 

0.0 

-4. 

0. 

n.O 

0.0 

-0. 

0. 

n.O 

0.0 

-0. 

0. 

i.O 

0.0 

-0. 

0. 

n.O 

0.0 

-0. 

ft. 

n.O 

0.0 

-0. 

4. 

n.O 

0.0 

-4. 

0. 

i.O 

O.n 

-0. 

0. 

n.O 

0.0 

-0, 

0. 

1.4 

4.0 

-4. 

0. 

1.0 

0.0 

-0. 

0. 

1.0 

O.o 

-4. 

0. 

n.O 

0.0 

-4, 

0. 

i.O 

0.0 

-4. 

0. 

n.O 

0.0 

-4. 

0. 

n.O 

0.4 

-4. 

0. 

1.0 

0.4 

-0. 

0. 

1.0 

0.4 

-0. 

0. 

1.0 

0.0 

-4. 

4. 

n.O 

0.4 

-4. 

0. 

n.O 

J.J 

-0. 

0. 

n.O 

0.4 

-0. 

0. 

n.O 

0.0 

-0. 

0. 

n.O 

4.0 

-0. 

0. 

n.O 

0.0 

•4. 

0. 

n.O 

4.1 

-0. 

0. 

n.O 

0.0 

-0. 

0. 

A.O 

0.0 

-0. 

0. 

n.O 

0.0 

-0. 

0. 

1.0 

0.0 

-0. 

0. 

1.0 

0.4 

0. 

.  0. 
4. 
0. 
0. 


3b. 

4. 

0. 

0. 
60. 
60. 
53. 
33. 
63. 
63. 

ii. 
45. 

140. 
144. 
115. 
114. 
73, 

25. 

0, 

4, 


.0. 

.0. 
.0. 


?47. 

1  n. 

1*1 . 

in. 

133. 

*. 

1*6. 

.9. 

1*5. 

•).. 

417. 

-JS. 

1.7. 

•411. 

.0. 

1). 

.4. 

1. 

.0. 

1. 

>S7. 

-411. 

U». 

-Si. 

174. 

•  rs. 

41?. 

-  '*. 

13. 
141. 

l?l  . 
191. 
1*4. 
474. 
**». 

*?S. 

4?7. 

4)9. 

.0. 
.1. 
.0. 


0. 

0. 
0. 
4. 
0. 
0. 
0. 
3. 
0. 


0. 
0. 
0. 

). 

0. 


.0. 
.0. 
.0. 
.0. 
.0. 
.0. 
.0. 
.0. 


41. 

*n. 


4«u    pOHrc4»l    Pnsir    (mn4 
kilt    iIGh(    A9!>LC    E4Q19 
kilt    imTcSifT    »4G9ITij'jF    ERRO* 
«<j    iNTmsili   ii»s 

9UH8tl*     Trf     '"39£C4S?S 


4.  l  riaf  :»srs 


34-99         48-9-1         73-1* 


?3. 

1  14. 

266. 

376. 

16. 

7b. 

140. 

1«9. 

3. 

1  4. 

22. 

14, 

-I. 

-9. 

-19. 

-  1  *  . 

11 

2* 

23 

?1 
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Table  4.  Best  track 


T'rt'vN  IRVING 


a. 5  r  tw»r« 


•HJ/J4/M. 

06o/ 1  2  ' 
oao  flu.' 

oao  so  V 
uS'.i  jti : 

J1J31  2.' 
u"jnn 
J  o  j  y  j  u  I 
jbu  y  it>, 

U91 Uou/ 
aal job/ 
091 j 12/ 
obi j i B/ 
uoi  100/ 
Obi  us; 

obi  1 12/ 

091  1  IB? 

OBl^OO/ 
081 206/ 

o b i  2 i  n 

0912197 
061 JU07 
UB1  306/ 

oai  jib/ 

joliui' 

-    -    ' .' 


"JilT 
I ». J     1  >7. 

I  3.0       I    18  . 


i 


■    !    ■ 


- 
-    >l 


-I    /  J   . 

J  4  i  /  J  e  I 

J  9  I    •  !       / 
jOI  '.8/ 

J8I SOU/ 
J  B  1  B  J  O  / 


1  3.. 
I  B.» 
IB.  J 
1B.0 
17.2 
IB.  3 
16.  4 

L7.S 

17. a 

19.4 
IB. 7 
19.2 

'0.0 

M.I 

'<r.j 
'^.  i 

'3.3 

-».o 

'6.2 


B.5 
-".3 


►  .  P 
.11.1 
1   '  J  .  P 

IH.l 
1  >9.  7 
1  '«.» 
l-'S.S 
1  p9.o 
1  '«.» 
1  '7.* 
I  '7.  I 
1  '6.  J 
13*,.* 
1  36.  7 
126.* 
l'6.n 
1'5.3 
I  'S.n 

:  3*.„ 
i  '».s 

IP*.. 


13. B 
I  •  .  I 
39.3 


:  'J.  > 

:  P3.  ' 

1?*. n 

:").fl 
">.  7 

1  ?9.g 

m.« 
i  u.s 

1  15.1 


0.0 

o.n 
o.n 
o.n 
o.n 
o.o 

1  '.a 

14.  1 
IB.  1 

lB.n 

15.  3 
19.8 
1  B.  7 
IS.  7 
1  7.3 
17. n 

17.? 

17.4 

17.5 
IT.  9 

18.  5 
14.  • 
20.1 
21.1 
22.0 

21.2 

23.  7 

24.1 

2*.* 

?6.3 

-3.  7 

<">.S 

I      ^ 

-B.  ) 

29.3 

31.1 

il.3 

32. a 

34.1 
Jb.3 
36. B 
39.3 
42.1 
44.  a 


O.n 

O.n 

o.n 

1  36. n 

I  36.1 

I  35.3 

I  35. n 
I  33. s 
I  32.7 
111.* 
I  30.  I 
128.  7 
I  27. « 

i  29.3 

1  ?B.« 
I  28.3 
127. S 
I  26. S 
I  25. » 
1  26.'" 
I  26.* 
1  36.  * 
1  ?S.7 
1  25.1 
I  ?*•« 
1  24. B 
I  ?».S 

l?».4 

12*.* 
I  '3.9 
i  23. « 
\23.7 
123. a 
123.7 
I  23.* 
I  ?*•» 
I  26.3 
I  29.* 
I  31  .* 
I  3*.l 
I  37. n 


30. 
Ju. 
30. 
30. 
25. 
25. 
30. 
30. 
30. 
30. 
35. 
*". 
55. 
35. 
35. 


'0. 

'0. 
'0. 

'0. 
75. 
90. 
85. 
^0. 

»o. 
■»o. 

■*u. 
''5. 
*0. 
90. 
BO. 
'0. 
60. 
30. 
3u. 
25. 
25. 


nsi 
-". 
-o. 
-o. 
-o. 
-o. 

-0. 

s. 
23. 
2  3. 
49. 
1  3. 
37. 
29. 
4*. 

8. 

50. 
21  . 
33. 
35. 
34. 
26. 
61. 

8. 
11. 
33. 
*!  . 
1  3. 

a. 
16. 

a. 
13. 
23. 
II. 
13. 

6. 
30. 

6. 
10, 

8. 
J*. 
33. 
36. 
36. 
«1. 


I|I9S 
4*  I  9  1 

0  . 
0. 
(1 . 
0. 

o. 

0. 
0. 

0. 


-6. 

0  . 

0. 

0. 
-5. 
-5. 
-5. 

5. 

0. 

0  . 

-5. 

0. 
0. 
0. 

o. 


5. 

0. 
0. 
5. 

0. 


►">'!  I 


O.n 

O.n 

O.n 

O.n 

O.n 

O.n 

19.3 

31  .n 

19.7 

I  9.« 

19.* 

PO.  I 

19.  ■ 

19.1 
1  7.1 
I  7.3 

18.3 

19. n 
18.7 
19. n 

19.4 

19. a 
33. o 
?5.n 
35.9 
37.* 
37. a 
37.. 
30. a 

PS. 3 
38.' 
>9.  i 
10. ~ 
31  .* 

13.' 

16.  1 

10.  < 

37.1 

18.7 

39.  I 

O.n 

O.n 

O.n 

O.n 


n ,  ii 

133.4 

13".  ■ 

I  3  <.  • 

1    34.4 

1  3n.s 

!?».- 
127.  . 

12s.  i 
12*.  * 
121.  ' 
<2s.- 
>24.3 

'2..- 
i2  >.. 

123.  - 

I2«..n 
I2n.» 
I27.n 
127.1 
12*. 4 
|2«..3 
12*.  t 

I  24,* 

|2*.« 

124.3 

12'.* 
I  2  <  .  ■• 

'  ?4  .  -» 

24..* 
I2*.s 


■  1  9r> 
0. 


0. 

). 
35. 


«ec»3i 

J3 '     •!• 


1  U. 
1  ho. 


"0. 

B5. 
'"0. 
1KB. 

I  uO. 

95. 
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MAPI.  East  China  Sea  Bottom  Contours 


SOUNDINGS    IN    FATHOMS   AND    FEE! 
;epth  contours  <h  ^rr-bMs 
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The  dispersive  arrival  at  the  prediction  site  of  the 
waves  propagating  from  a  single  source  region  is  repre- 
sented by  a  plot  of  period  versus  arrival  time  in  Figures 
9,  10  and  11.   The  energy  in  each  low  frequency  (7~19  sec) 
spectral  component  of  the  swell  at  Cheju-do  is  calculated 
by  modifying  the  energy  spectrum  at  the  source  for  the 
refraction  and  shoaling.   The  lowest  period  (T  =  4  sec) 
is  assumed  totally  eliminated  by  attenuation. 

The  total  propagated  energy  at  any  given  time  is 
estimated  by  summing  all  the  directional  components. 
The  predicted  swell  heights  H 1/10  are  assumed  to  be  related 
to  the  energy  by  H  1/1 0  *k /E ,  where  K  is  constant  (0.4)  and 
E  is  energy  in  a  unit  of  2  •  kiloerg/cm   and  H  1/10  is  in 
meters . 

By  graphing  the  energy  associated  with  the  5   mean 
period  bands  as  a  function  of  arrival  time  at  the  forecast 
site  and  summing,  then  a  plot  of  the  total  propagated 
energy  is  obtained  (see  Figures  13,  15  and  17).   This 
plot  is  compared  with  the  observed  wave  heights  in 
Chapter  IV-D. 

E.   TERMINOLOGY 

To  make  clear  the  notation  used  in  this  study,  the 
following  terms  and  symbols  are  defined. 
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sea   -  Those  waves  found  within  a  generating  area, 
swell  -  Those  waves  which  have  moved  out  of  a 

generating  area. 
RQ    -  The  distance  from  the  source  to  the  "bottom 

line"  where  the  corresponding  periods  feel 

"bottom"  i.e.,  the  swell  travel  distance 

in  deep  water.   The  depth  at  the  "bottom 

line"  is  h  =  !^. 
R     -  The  distance  from  the  "bottom  line"  to  the 

forecast  site. 
T     -  The  mean  period  in  a  spectral  band  (T  =  1/f) 
Hl/10   -  The  average  heights  of  the  1/10  highest 

waves. 
H  1/3  -  The  average  heights  of  the  1/3  highest 

waves  (significant  wave  height) . 
tar   -  The  swell  arrival  time  (GMT) . 
t     -  The  swell  travel  time. 
L     -  Wavelength  (L   in  deep  water) . 
Cg    -  Group  velocity  (Cg   in  deep  water) . 
C     -  Phase  speed  (Co  in  deep  water) . 
D     -  The  total  swell  travel  distance  (D  =  RQ  +  R) 
g     -  Acceleration  of  gravity 
g     -  Wave  angular  frequency  (2/C/L)  . 
k     -  Wave  number  (27C/L)  . 
h(d)   -  The    ater  depth. 
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ao     "  The  an^le  between  wave  crest  and  bottom  contour 

at  the  depth,  h  =  5°  . 

2 

T4     -  The  total  energy  at  grid  point  4  in  the  generat- 
ing area. 

D4     -  The  directional  energy  component  at  grid  point  4 
in  generating  area. 

Wave  spectrum  (energy  spectrum)  -  The  distribution  of 
either  wave  height  or  energy  with  period.   The 
potential  energy  of  the  sea  surface  is  proportional 
to  the  mean  of  the  square  of  the  wave  height. 

Dispersion  -  A  process  which  leads  to  longitudinal 

spreading  of  the  wave  energy  as,  in  deep  water, 
energy  in  each  spectral  component  is  propagated 
at  a  characteristic  group  velocity  Cg,  the  long 
waves  having  the  larger  group  velocity. 
-  The  mean  directional  bands  of  16  unit  points 
rose. 

H'     -  The  wave  height  before  refraction  in  deep 
water. 


24 


II.   PREDICTION  OF  TROPICAL  STORM  WAVES  BY  TYWAVES 

A.   LOCATION  OF  THE  SELECTED  POINT-SOURCES 

The  TYWAVES  program  produces  a  printout  of  the  complete 
energy  spectrum  for  12  points  around  each  typhoon.   Each 
point  lies  at  a  chosen  distance  from  the  storm  center  in 
a  fixed  direction,  no  matter  what  the  storm  movement  direc- 
tion is.   These  points  are  representative  origins  of  the 
wave  energy  emerging  from  the  typhoon. 

Only  those  points  are  selected  which  are  in  the  region 
of  maximum  winds  and  where  important  energy  components  are 
directed  toward  the  distant  prediction  site. 

The  arrangement  of  points  whose  wave  spectra  are  used 
in  this  study  is  shown  below. 

+  9 

+5     +1     +6 

+11   +3     t#    +4   +12 

+7     +2     +8 
+10 

The  distances  from  the  storm  center  to  points  1,  2,  3 
and  4  are  three  grid-lengths  (3  x  40  =  120  NM) ,  to  9,  10, 
11  and  12  are  five  grid  lengths  (5  x  40  =  200  NM) ,  and  to 
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5,6,7  and  8  are  3  x  /2  grid  lengths  (3  x  /2  x  4  0  ---   17  0  NM)  , 
respectively. 

Among  the  12  points,  only  points  4,6,8  and  12  were 
selected  as  sources,  because  the  spectral  energy  from  only 
these  4  points  can  propagate  to  Cheju-do,  Korea  (33.2°N 
126. 6°E) . 

The  locations  of  the  sources  for  each  time  of  interest 
for  the  three  storms  are  shown  below  in  Table  VI,  which 
was  derived  and  calculated  from  the  postanalysis  "best 
track"  data  shown  in  Tables  III,  IV,  and  V  [2].   The 
relationship  between  latitude  and  longitude  is: 

Lat.  =  Long,  x  cos  (Lat.) 

=  Long,  x  cos  { (storm  center  +  forecast)/2} 

lat 
e.g.: 

When  the  storm  center  (best  track)  is  20.0  N  126.7  E, 

point  4  is  same  lat  (20.0°N),  but  120  NM  (2°  in  lat) 

east  of  the  center. 

„o  T      oT           ,20.0  +  20.0, 
Thus,  2   Lat  =   Long  x  cos  ( ^ > 

O,  -I  A° 

=   Long  x  cos  10 

.*.  °  Long  =  2°  Lat/cos  10°  =  2.128°  E 

=  126.7  +  2.128 


=  128.8°  E 


Therefore,  point  4  *  20.0°N  128. 8°E. 
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B.   SEA-SPECTRA  AT  THE  POINT  SOURCES 

The  description  of  the  wave  fields  using  the  spectral 
component  method  developed  by  NEPRF  was  used  in  this  study. 
Only  those  directional  components  able  to  propagate  to  the 
forecast  site  need  be  considered. 

A  schematic  example  of  the  TYWAVES  spectral  prediction 
in  a  typhoon  area  is  given  below.   The  whole  spectral  table 
for  the  selected  sources  is  shown  in  Appendix  B. 
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T  -  6  spectrum 
(80  components) 


Total 
energy 


|T  -  spectrum 
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^period  bands 


Table  VII  is  one  example  of  a  period-direction  spectrum 
at  source  point  4  in  Typhoon  Owen,  at  06  on  2  5  September 
1979. 

The  mean  period  is  13  sec  with  the  maximum  energy  86 
(26-  kiloeng/cm  )  in  the  dominant  period.   The  dominant 
wave  direction  is  SE,  but  the  energy  from  that  direction 
will  not  arrive  at  the  forecast  site,  only  the  SSE  component 
will  arrive  at  the  site  (as  seen  in  Figure  5) . 
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C.   SELECTION  OF  THE  DIRECTIONAL  ENERGY  AT  POINT  SOURCES 

A  constant  energy  was  assumed  between  parallel  rays 
separated  by  width  80-160  NM  (2-4  grid  length),  which  is 
40^80  NM  right  at  the  grid  point  depending  upon  the  typhoon 
size.   This  is  why  the  energy  in  each  component  is  assumed 
constant  along  the  deep  water  propagation  path.   As  example 
of  possible  refraction  diagrams  under  the  above  considera- 
tion and  bottom  contours  in  Map  1  is  shown  in  Figure  5. 

Only  the  SSE  component  was  considered  to  propagate  to 
the  forecast  site.   The  SE  component  would  refract  toward 
China,  west  of  Cheju-do,  and  the  S  component  would  pass  to 
the  east  of  Cheju-do  as  seen  in  Figure  5. 

All  the  directional  components  were  chosen  in  the  same 
way,  and  are  underlined  as  shown  in  Table  7. 


>^>  si 

U        LINE  l\Vft\       J 

1 


ea: 


Figure  5.   Possible 
refraction  diagram  at 
point  4  at  06  GMT 
September  25,  1979. 
The  curves  around 
typhoon  center  are 
contours  of  wave  height 
(H  1/3  in  feet) .   The 
detailed  bottom  contours 
are  shown  on  Map  1. 
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TABLE  7 

Period-Directional  Spectrum  at  Point  4  at  06  GMT  September 
25,  1979.   The  energy  units  are  26-  kiloerg/cm2. 

TYPHOON    WAVES       3C I  IT  U. 
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D.   TOTAL  AND  DIRECTIONAL  ENERGIES  AT  POINT  SOURCES 

Figures  6 ,  7  and  8  show  the  energy  components  at  each 
selected  grid  point  in  the  typhoon  area.   Each  line  is 
labeled  by  selected  grid  point  (source  point) .   T   indicates 
total  energy  given  by  TYWAVES  at  point  n  and  D   represents 
the  sum  over  all  periods  of  the  energies  with  proper 
direction  to  reach  Cheju-do. . 
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V9o>    =il1  Ei    <V    V 

For   example,    where   n   =    12,    0    is   SE,    then 

D12(9=SE)     =    E?     (T=4,     9q=SE)     +    E    Q     (T=10,     9q=SE) 

+    E13     (T=13,     6q=SE)     +    E16     (T=16,     eQ=SE) 

+    E19     (T=19,     9    -SE) . 

This  represents  the  total  energy  at  point  12  directed 
from  SE  toward  Cheju-do. 

The  energy  is  in  practice  converted  into  height  H  1/10 
by  making  H  1/10  =  0.4  /I,  which  is  described  in  Chapter  I.D. 
Some  special  features  of  the  wave  fields  in  each  source  are 
now  discussed. 

1.  Typhoon  Hope 

The  maximum  sustained  wind  speed  was  over  100  kts 
for  two  days  from  00  GMT  July  31.  The  maximum  wave  height  (H  1/10) 
recorded  was  12.80  m.   The  directional  energy  components  from 
grid  point  12  were  ignored  from  00  GMT  July  31  to  00  GMT 
August  01,  because  they  were  too  small. 

2 .  Typhoon  Irving 

Only  the  propagation  of  energy  components  until 
00  GMT  August  16  was  considered,  because  the  forecast  site 
is  inside  of  typhoon  area  after  that  time.   Therefore,  after 
00  GMT  August  16,  the  propagated  wave  energy  should  be 
combined  with  local  energy  at  prediction  site. 
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Among  the  directional  components,  only  those  from 
grid  points  4  and  12  were  propagated  to  the  prediction  site 
during  the  typhoon  period,  because  the  energies  from  other 
sources  were  refracted  away  from  Cheju-do. 

3.   Typhoon  Owen 

After  06  GMT  September  29,  the  typhoon  passed  east 
of  southern  Japan.   No  further  propagation  to  Cheju-do 
occurred. 

In  most  cases,  only  one  directional  energy  component  was 
selected  at  each  grid  point.   But  both  components  (SE  and 
SSE)  were  able  to  reach  the  forecast  site  on  18  GMT  September 
27,  from  grid  point  4,  because  of  the  width  of  their  ray 
boundaries  at  distances  near  the  source. 
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III.   PREDICTION  OF  SWELL  CHARACTERISTICS 

The  periods  of  high  waves  due  to  typhoons  are  usually 
in  the  range  of  7~19  sec.   They  travel  in  deep  water  with 
the  group  corresponding  velocities  10-30  kts.   It  normally 
necessitates  1-4  days  for  conspicuous  swells  generated  in 
a  typhoon  in  the  far  south  seas  to  reach  Cheju-do,  Korea. 

The  group  velocities  of  periods  7  and  10  sec  are 
constant  along  their  paths  since  they  remain  in  "deep 
water."   But  for  periods  13,  16  and  19  sec,  the  group 
velocities  vary  because  of  depth  variations.   Thus,  I 
considered  two  distance  components:   in  deep  water  R  and 
in  shallow  water  R  for  periods  13,  16  and  19  sec. 

The  travel  distances  between  grid  points  and  observa- 
tion points  are  calculated  with  the  assumption  of  plane 
geometry.   Then,  group  velocities,  travel  distances  and 
travel  times  are  derived  for  each  period  in  this  section. 

A.   THE  SETTING  OF  THE  FORECAST  SITE  AND  ITS  RELATION  TO 
THE  WAVE  SOURCES 

The  passes  of  the  Ryukyu  Islands  act  as  windows  between 
the  energy  sources  and  the  forecast  site,  Cheju-do,  Korea. 
Since  the  largest  land  length  for  blockage  of  energy  is 
only  about  40  NM,  this  length  is  not  sufficient  to  inter- 
rupt totally  the  energy  propagation.   Thus,  as  a  further 
simplification,  I  neglected  the  effect  of  the  Ryukyu  Islands 
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against  the  energy  propagation.   The  window  is  assumed 
sufficient  for  total  propagation. 

B.   GROUP  VELOCITY 

Some  simplications  are  used  in  assigning  group  velocity 

values  along  the  route  from  typhoons  to  the  observation 

site. 

From  linear  wave  theory  group  velocity,  C   is  given  by 

y 

Cg  "  2  [C  +  sinh  2kh]   =  nc 

(all  the  symbols  are  defined  below) 
and  phase  speed  c  is  given  by  c  =  -  tanh  kh.   In  deep  water 
this  is  well  approximated  by  c   =  g/a  =  3.03  T  (ktc)  for 
T  given  in  seconds,  and  in  shallow  by  c  =  /gh.   But  in  the 
general  case  c  varies  with  both  the  depth  of  water  and  the 
wavelength. 

The  classifications  "deep"  and  "shallow"  are  given  in 
terms  of  "relative  depth  h/L  ,  described  below.   Since 
there  is  no  "shallow  water"  between  the  wave  sources  and 
Cheju-do,  Korea,  (corresponding  to  the  periods  used  by 
TYWAVES) ,  I  will  use  only  the  general  equation  and  the 
deep  water  approximation.   In  these  equations 

C  =  group  velocity 

c   =  phase  speed  (c   in  deep  water) 

n   =  the  fraction  of  energy  propagated  at  phase 
speed 
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g   =  gravity  acceleration 

h  =  water  depth 

a   =  wave  angular  frequency  (2T/T  =  2kf,  where  T  is 

period) 
k  =  local  wave  number  (2W/L,  where  L  is  wave  length) 
f   =  wave  frequency  (1/T) 
It  is  important  to  note  that  c  and  group  velocity,  Cg 
must  be  found  as  functions  of  both  h  and  a.   A  straight- 
forward method  is  illustrated  below,  where  the  dependence  on 
kh  is  replaced  by  the  deep  water  wavelength  L  =  2T~  g/T  . 

The  local  wavelength  is  then  found  from  L/L  ,  a  function 
of  the  relative  depth  h/L  . 

Since  C/C   =  L/L  ,  the  phase  speed  c  can  be  found  from 
h  and  a. 

To  estimate  the  group  velocity  by  definition,  C   = 

y 

i  c  [1  +  — « — r-nr]   =  nc  where,  n  takes  on  the  following 
2    L     smh  khJ 

values  for  corresponding  h/L  : 

general  case  (y  >  h/L   _>  jq)  :  l>n>y  , 

deep  water  case  {h/L     >  y) :  n  =  y  • 

In  order  to  simplify  calculations,  I  used  n  =  3/4  in 
the  waters  between  the   Ryukyu  Islands  and  Cheju-do 
where,  for  the  longer  periods,  the  relative  depth  range 

1  >h/L  >-W  .   With  this  simplication,  the  corresponding 

2  ~    o  20 

group  velocities  become  the  following: 
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T  =  13  sec;  C   =  2  9.4  kts 

y 

T  -  16  sec;  C   =  33.6  kts 
T  =  19  sec;  C   =  35.8  kts 
The  influence  on  the  calculation  of  C   of  the  use  of 

g 

these  representative  constant  values  is  discussed  in 
Chapter  IV. D. 

For  deep  water 

Cq  =  3.03T  (kts  for  T  in  sec) 

C  =  nc  =  T  C   =  1.515T  (kts  for  T  in  sec) 
y  o 

Therefore,  group  velicity  in  knots  and  depth  at  which 

C   replaces  C   are  shown  below: 

T  (sec)  7      10       13  16        19 

C   /C   (kts)   10.6/-   15.2/-   18.7/29.6   24.3/33.6   28.8/35.C 
go   g 

Lo 
h(-j-  in  ft)    126      256        433  656       924 

C.   TRAVEL  DISTANCE  AND  ARRIVAL  TIME 

The  travel  distance  between  grid  points  in  the  storm  and 
the  observation  point  33.2  N  126.6  E  are  calculated  with  the 
assumptions  of  plane  geometry.   Consider  this  example. 

On  12  GMT  July  30,  the  grid  point  4  was  at  17. 4°Nl33 . 9°E. 

Thus,  the  north-south  component  is  given  by 
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y  =  [lat  (site)c  -  lat  (grid  point)0]  x  60  MM 
=  [33.2  -  17.4]   x  60 
=   94  8  NM 
and  the  east-west  component  is  approximated  by 

x  =  [long  (grid  point)  -  long  (site)]   cos  (mean  lat)] 

x  60  NM 
=  [133.9  -  126.6]   cos  (17.4)  x  60  NM 
=  418  NM 

Thus,  the  travel  distance,  D,  is  given  by 

2     2  ^ 
D  =  (x   +  y  )  2  and  for  grid  point  4, 

D  =  (4182  +  9482)^  =  1036  NM 

The  travel  time,  t,  is  given  by 

R  R 

t  =  _£  +  7r-   ,  where  R  +  R  =  D 

c  c  o 

9      g 

^o     3 

R  =  the  deep  water  distance 

R  =  the  remaining  distance 

Therefore,  the  arrival  time,  tar  of  each  period  band  is 

given  by 

tar  (T)  =  t   +  t,  where  t   =  leaving  time  from 
o  o         3 

typhoon  area 
t   =  travel  time 
The  Figures  9,  10  and  11  show  the  swell  arrival  time 
for  each  of  the  three  typhoons.   Each  curve  indicates  the 
arrival  time  of  energy  from  a  source  at  a  specified  leaving 
time  labeled  with  various  periods  (from  7  to  19  sec). 

40 


Appendix  C  shows  the  arrival  time  calculations  for 
each  typhoon. 
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IV.   THE  ENERGY  PROPAGATION  PROCEDURE 

A.   PROPAGATION  OF  SPECTRAL  ENERGY  COMPONENTS 

The  spectral  forecast  permits  the  forecaster  to  deal 
with  only  that  range  or  periods  which  have  important  energy. 
Each  spectral  component  is  tracked  with  its  respective  group 
velocity,  only  those  directions  being  chosen  for  which  waves 
can  reach  Cheju-do.   When  the  typhoon  is  moving  with  a 
velocity  component  toward  Cheju-do,  only  spectral  components 
with  a  group  velocity  greater  than  the  movement  of  the 
typhoon  are  considered  (as  in  Typhoon  Irving) . 

Wave  energy  generated  in  a  relatively  small  region  at 
all  frequencies  will  spread  over  a  much  larger  region  as  it 
propagates  outward  from  its  source,  and  the  wave  character- 
istics change  in  such  a  way  as  to  become  more  "swell  like." 
There  are  essentially  three  processes  which  contribute  to 
this  change  in  wave  characteristics:   dispersion  and  angular 
spreading,  which  are  modified  by  shoaling  and  refraction, 
and  attenuation. 

In  the  procedure  followed  here,  I  have  ignored  attenua- 
tion for  waves  of  long  period.   Dispersion  and  angular 
spreading  are  accounted  for  by  simply  following  components 
to  Cheju-do  with  appropriate  shoaling  and  refraction 
factors. 
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For  ease  of  calculation  shoaling  and  refraction  processes 
are  simplified  as  described  below.   As  seen  on  Map  1,  a 
slightly  curved  contour  of  100  fathoms  (590  ft)  connects 
the  northern  tip  of  Taiwan  to  the  south-western  tip  of 
Japan.   The  slope  from  590  ft  bottom  line  to  413  ft  (70 
fathoms)  is  very  steep.   But  most  bottom  topography  along 
the  path  to  Cheju-d0/  north  of  the  413  ft  bottom  contour 
is  almost  flat  with  the  depth  of  45  fathoms  (266  ft)  up  to 
the  forecast  site  (33.2°N  126. 6°E). 

Therefore,  I  have  approximated  to  underwater  topography 
by  assuming  only  two  water  depths,  a  deep  water  and  a 
shallow  water  (intermediate  water  depth)  region  with  an 
abrupt  jump  between  them. 

To  compute  the  wave  characteristics  at  a  shoal  water 
site,  the  shoaling  and  refraction  are  considered  using 
the  values  of  Cg  and  n  from  Chapter  III.B. 

The  energy  of  component  Eq  (t,9)  in  the  typhoon  area 
(deep  water)  is  transformed  after  shoaling  and  refraction 
to  its  value  at  the  forecast  site  according  to 
Ed,  J)  =  Eq  (i,  j)-k2  (±#  -j  k2  {if    j} 
where  E  (i,  j)  is  the  energy  of  the  component  of 
period  T\  and  which  had  the  direction  9 .  in  the 

-J  • 

generating  area, 

and  ks  (i,  j),  kr  (i,  j)  are  the  respective 

shoaling  and  refraction  factors  of  those  components. 
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Details  of  calculation  of  k   and  k   are  given  in  the 

S         r 

following  sections.   The  refraction  and  shoaling  calcula- 
tions for  three  typhoons  are  shown  in  Appendix  D.   The 
following  is  a  sample  calculation  on  12  GMT  July  3  0  for 
Typhoon  Hope  at  source  point  4. 

T         7         10         13         16         19 
ks        1  1         0.92      0.89       0.90 

kr  1  1                     1  0.989      0.972 

(ks  kr)  1  1         0.8464  0.7748     0.7653 

EQ  7  19  11  9  14 

E  7  19         9.3  7.0  10.7 

where  k   is  derived  from  plate  C-l  [6] 

with  h  =  40  fathoms  =  2  36  ft 

k   is  derived  from  Figure  2-19  [5] 

with  h  =  413  ft  and  a     =   40°. 

o 

B.   SHOALING  AND  REFRACTION  FACTOR 

1.   Shoaling  Factor 

The  shoaling  factor,  k  =  H/H '  was  derived  from  the 
3  so 

plate  C-l  [6]  with  depth  and  period  of  the  appropriate 
spectral  component.   The  shoal  site  depth  near  Cheju-do 
(33.2°N  126. 6°E)  is  40  fathoms,  and  the  relative  depth  h/L 
and  shoaling  factor  of  each  period  are  shown  below. 


47 


10 
8.0 


T-rq 


0.00 


0.0001  0.0002  0.0004     0  00. 0.002   0.004     0.01    0.02    0.04  0.06   0.1     0.2 


0.4  0.6    1.0 


d 
Lo 


iftM  Wiegel.  R.L.,  "OitilUtory  W.»c"  U.S.  Amy,  fl««ch  Eroaon  Bo^rd. 
DuJcon,  Special  Imu*  No.  1.  Juiy  1948. 


Plate  C-l     Illustration  of  Various  Functions  of  -£- 


48 


V 

Q. 
O 

to 

e 
o 

c 
o 


o: 


s 

3 

Q 

—    i/> 

1>    i- 

X§ 
T3    — 

ss 

=  u 

2  * 

=  Q 
Q  — 

»   — 

o  o 
5  3 

c  Q. 

»/>  *- 

4)  -c 

Ss 

-C    — 


0> 


3 
0< 


49 


10 

13 

16 

19 

0.46 

0.24 

0.18 

0.13 

1 

0.92 

0.89 

0.90 

T       17 
h/LQ    0.94 

ks      1 

2 .   Refraction  Factor 

Generally,  two  basic  techniques  of  refraction 
analysis  are  available — graphical  and  numerical.   Several 
graphical  procedures  are  available,  but  all  methods  of 
refraction  analyses  are  based  upon  Snell's  law.   Refraction 
may  be  treated  analytically  in  any  region  with  straight 

and  parallel  contours,  by  using  Snell's  law  directly: 

c 
sin  a =  (— )  sin  a  ,  where  a  is  the  angle  between  the  wave 

co        ° 
crest  and  the  bottom  contour,  and  aQ  is  the  angle  between 

the  deep  water  wave  crest  and  the  bottom  contour. 

Figure  2-19  [5]  shows  the  relationships  among 

a,  a  ,  period,  depth  and  refraction  factor  in  graphical 

form.   I  derived  the  refraction  factor  from  using  this 

graph,  Figure  2-19,  the  bottom  contour  and  period.   Also, 

I  assumed  that  the  refraction  factor  is  1  if  the  angle  a 

3     o 

(between  crest  and  bottom  contour)  is  less  than  10   and 
that  the  refraction  occurs  only  one  time  at  the  depth  of 
maximum  bottom  gradient  where  waves  are  refracted  because 
of  the  wide  flat-bottomed  portions  over  most  of  the  inter- 
mediate water  propagation  path  as  seen  on  Map  1. 
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With  these  considerations,  I  derived  the  refraction 
factor  and  predicted  the  shoaling  energy.  Those  procedures 
are  shown  in  Appendix  D. 

C.   SHOALING  AND  REFRACTION  OF  THE  SPECTRAL  COMPONENTS 

Each  deep  water  wave  spectral  component  derived  from 
the  TYWAVES  model  moves  with  its  respective  group  velocity. 
I  considered  the  80  energy  components  to  behave  as  mono- 
chromatic component  waves.   To  assess  the  shoaling  and 
refraction  of  each  to  the  shoal-water  site  I  assumed  that 
the  wave  power  transmitted  between  a  given  pair  of 
orthogonals  remains  constant  at  all  depths  (this  means  no 
frictional  losses,  diffraction  or  scattering,  and  also 
implies  that  a  steady  state  exists) . 

With  these  assumptions  the  wave  power  P  is  given  by 
P  =  EC   b=E   C    b 

g     o  gQ  o 

C    bQ 
Thus,  E  =  Eq  _!o_^  =  Eq  x  k2  k2  m    1  pgH2 

g 

where  E  =  the  average  wave  energy  per  unit  area 

of  sea  surface   for  waves  transformed  by  shoaling  and 

refraction. 

P  =  water  density 

g  =  acceleration  of  gravity 

H  =  wave  height  of  transformed  waves 

C   C    group  velocity 
gQ'  g 
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bQ,b  =  orthogonal  separations 
k   =  shoaling  factor 
kr  =  refraction  factor 

The  energy  in  each  spectral  component  of  the  swell  at 

the  observation  site  was  calculated  by  modifying  the  energy 

spectrum  at  the  selected  sources  for  the  effects  of  shoaling 

and  refraction  according  to  equation  E  =  E  -k2-k  . 

o   s   r 

The  energy  associated  with  the  various  components  as  a 
function  of  time  of  arrival  at  observation  site,  as  seen  in 
Figures  9,  10  and  11. 

The  total  energy  in  the  swell  at  any  given  arrival  time 
is  estimated  by  summing  all  the  shoaling  components  at  that 
time. 

In  summary,  each  component  in  typhoon  area's  T  -  6 

spectrum  is  shoaled  and  refracted  using  the  k   and  k 

3       s       r 

values  appropriate  to  it  to  find  the  energy  at  the  observa- 
tion site. 

1.  Shoaling  Process  (Computation) 
See  Appendix  D. 

2 .  Shoaling  Energy  Components  From  Each  Source  Versus 
Arrival  Time 

As  seen  in  Figures  12,  14  and  16,  the  energy  of 

the  components  from  all  sources  is  shown  as  a  function  of 

its  arrival  time  at  Cheju.   See  Appendix  D. 
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3.   The  Predicted  Swell  Waves  at  Prediction  Site 

As  shown  in  Figures  13,  15  and  17,  the  predicted 
swell  waves  are  the  sum  of  all  transformed  components  at 
given  time  at  the  forecast  site. 

For  Typhoon  Irving  (see  Figure  15)  beginning 
12  GMT  August  16,  the  forecast  site  is  already  inside 
typhoon  area.   Therefore,  there  is  no  prediction  done 
after  that  time.   These  predictions  are  discussed  in 
Chapter  IV. D. 

D.   THE  OBSERVED  DATA 

Observations  of  wave  conditions  for  verification  of 
the  swell  forecasts  were  obtained  from  the  sources,  the 
National  Oceanographic  Data  Center  (NODC)  and  the  Republic 
of  Korea  (ROK)  Navy.   All  listed  heights  are  based  on 
visual  estimates. 

Following  Table  VIII  shows  the  visual  observation  of 
swell  heights  (H  1/10)  for  each  typhoon.   The  data  from 
NODC  are  sparse  and  often  far  from  the  forecast  site, 
but  they  appear  to  be  samples  from  the  same  set  as  those 
of  the  ROK  Navy.   The  data  from  ROK  Navy  visual  observa- 
tions were  made  at  33.2  N,  126.5  E  close  to  my  point  of 
interest  (33.2°N,  126. 6°E).   So  I  considered  this  point 
exactly  the  same  as  my  forecast  site. 
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TABLE  8 
The  Visual  Observation  Data  by  NODC  and  ROK  Navy 
For  Typhoon  Hope 


DTG  ( GMT ) 

Location  (  N-  E) 

H  1/10  (meters) 

Source 

79080100 

33.6-125.2 

2.5 

NODC 

i 

0104 

32.3-126.0 

2.5 

NODC 

I 

0112 

31.2-126.7 

4.0 

NODC 

I 

0118 

27.7-124.6 

6.5 

NODC 

i 

0121 

27.8-125.4 

6.5 

NODC 

I 

0200 

29.1-128.0 

6.5 

NODC 

I 

0206 

28.0-129.0 

5.5 

NODC 

i 

0218 

26.9-125.8 

5.0 

NODC 

i 

0300 

28.7-128.6 

3.0 

NODC 

i 

0306 

31.4-126.2 

3.0 

NODC 

I 

0312 

32.3-125.9 

3.5 

NODC 

I 

080100 

33.2-126.5 

3.0 

ROK 

Navy 

0104 

33.2-126.5 

3.5 

ROK 

Navy 

0108 

33.2-126.5 

3.7 

ROK 

Navy 

0112 

33.2-126.5 

4.0 

ROK 

Navy 

0116 

33.2-126.5 

4.5 

ROK 

Navy 

0120 

33.2-126.5 

4.5 

ROK 

Navy 

0208 

33.2-126.5 

5.0 

ROK 

Navy 

0212 

33.2-126.5 

4.5 

ROK 

Navy 

0216 

33.2-126.5 

5.5 

ROK 

Navy 

0220 

33.2-126.5 

4.5 

ROK 

Navy 

0308 

33.2-126.5 

3.5 

ROK 

Navy 

0312 

33.2-126.5 

3.5 

ROK 

Navy 

0316 

33.2-126.5 

3.0 

ROK 

Navy 

0320 

33.2-126.5 

5.0 

ROK 

Navy 

0408 

33.2-126.5 

4.0 

ROK 

Navy 

0412 

33.2-126.5 

3.0 

ROK 

Navy 

0416 

33.2-126.5 

3.0 

ROK 

Navy 

0420 

33.2-126.5 

2.5 

ROK 

Navy 

60 


2.   For  Typhoon 

Irving 

DTG  (GMT) 

Location  (°N-°E) 

H  1/10  (meters) 

Source 

79081400 

29.10127.6 

2.5 

NODC 

1403 

28.5-127.6 

4.0 

NODC 

1406 

28.7-129.3 

3.5 

NODC 

1409 

27.7-129.7 

4.0 

NODC 

1412 

29.2-131.1 

4.5 

NODC 

1500 

30.0-125.0 

6.0 

NODC 

1503 

30.7-127.3 

4.0 

NODC 

1515 

28.3-130.3 

8.0 

NODC 

1518 

27.2-130.0 

9.0 

NODC 

1521 

26.5-128.8 

8.0 

NODC 

1600 

27.1-130.0 

6.5 

NODC 

1603 

27.3-128.5 

8.0 

NODC 

1606 

29.6-129.2 

5.5 

NODC 

1700 

27.4-125.5 

6.0 

NODC 

I 

79081400 

33.2-126.5 

2.5 

ROK 

Navy 

1408 

33.2-126.5 

2.5 

ROK 

Navy 

1412 

33.2-126.5 

3.0 

ROK 

Navy 

1416 

33.2-126.5 

3.5 

ROK 

Navy 

1420 

33.2-126.5 

3.5 

ROK 

Navy 

1508 

33.2-126.5 

4.0 

ROK 

Navy 

1512 

33.2-126.5 

4.5 

ROK 

Navy 

1516 

33.2-126.5 

5.0 

ROK 

Navy 

1520 

33.2-126.5 

6.0 

ROK 

Navy 

1608 

33.2-126.5 

8.0 

ROK 

Navy 

1612 

33.2-126.5 

9.0 

ROK 

Navy 

1616 

33.2-126.5 

8.0 

ROK 

Navy 

1620 

33.2-126.5 

8.0 

ROK 

Navy 

1708 

33.2-126.5 

8.0 

ROK 

Navy 

1712 

33.2-126.5 

9.0 

ROK 

Navy 

1716 

33.2-126.5 

8.0 

ROK 

Navy 

1720 

33.2-126.5 

7.0 

ROK 

Navy 
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3 .      For  Typhoon   Owen 


DTG  (GMT) 

Location  (°N-°E) 

H  1/10  (meters) 

Source 

79092600 

31.4-128.0 

4.0 

NODC 

2606 

33.6-129.1 

3.5 

NODC 

2615 

31.8-129.5 

3.0 

NODC 

2700 

28.1-126.9 

6.0 

NODC 

2703 

27.1-126.6 

6.5 

NODC 

2706 

32.6-125.7 

4.0 

NODC 

2709 

31.6-125.7 

5.5 

NODC 

2712 

30.7-129.8 

6.5 

NODC 

2715 

28.7-125.3 

6.0 

NODC 

2800 

32.8-128.3 

5.5 

NODC 

2803 

30.6-125.7 

6.5 

NODC 

2806 

31.5-126.8 

6.5 

NODC 

2809 

31.6-126.6 

6.0 

NODC 

2812 

31.7-126.6 

6.5 

NODC 

2818 

34.0-128.3 

6.0 

NODC 

2900 

32.7-127.7 

6.5 

NODC 

2606 

32.1-127.3 

5.5 

NODC 

2912 

31.0-126.4 

5.5 

NODC 

2918 

34.0-129.7 

5.0 

NODC 

3000 

29.0-124.7 

4.0 

NODC 

3006 

30.8-127.7 

3.0 

NODC 

3018 

30.4-131.6 

4.0 

NODC 

79092600 

33.2-126.5 

2.5 

ROK 

Navy 

2608 

33.2-126.5 

4.0 

ROK 

Navy 

2612 

33.2-126.5 

5.0 

ROK 

Navy 

2616 

33.2-126.5 

5.0 

ROK 

Navy 

2620 

33.2-126.5 

4.5 

ROK 

Navy 

2708 

33.2-126.5 

5.5 

ROK 

Navy 

2712 

33.2-126.5 

5.0 

ROK 

Navy 

2716 

33.2-126.5 

5.0 

ROK 

Navy 

2720 

33.2-126.5 

5.5 

ROK 

Navy 

2808 

33.2-126.5 

4.5 

ROK 

Navy 

2812 

33.2-126.5 

5.0 

ROK 

Navy 

2816 

33.2-126.5 

5.5 

ROK 

Navy 

2820 

33.2-126.5 

5.5 

ROK 

Navy 

2908 

33.2-126.5 

4.5 

ROK 

Navy 

2912 

33.2-126.5 

4.0 

ROK 

Navy 

2916 

33.2-126.5 

4.5 

ROK 

Navy 

2920 

33.2-126.5 

4.5 

ROK 

Navy 

3008 

33.2-126.5 

5.0 

ROK 

Navy 

3012 

33.2-126.5 

4.0 

ROK 

Navy 

3016 

33.2-126.5 

2.5 

ROK 

Navy 

3020 

33.2-126.5 

3.0 

ROK 

Navy 
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E.   COMPARISON  OF  SWELL  PREDICTIONS  AND  OBSERVATIONS 

The  predicted  and  observed  heights  are  plotted  in 
Figures  13,  15  and  17,  and  are  good  agreement.   Detailed 
comparisons  for  each  typhoon  follow: 

1.  Typhoon  Hope 

The  first  predicted  swell  arrival  times  from  each 
of  three  sources,  and  the  NODC  and  ROK  Navy  observations 
are  almost  the  same  with  similar  energies  (see  Figure  13) . 
The  peak  energy  arrival  times  are  also  nearly  the 
same,  but  the  peak  predicted  energy  is  slightly  higher  than 
the  observations.   This  is  reasonable  agreement  since  it  is 
estimated  that  the  visual  observations  can  be  error.   Verploegh 

1961  [18]  estimated  the  average  observational  error  for  a 
visual  observation  of  wave  height  varies  from  1  ft  at  5  ft 
wave  heights  to  3  ft  at  18  ft  wave  heights. 

2.  Typhoon  Irving 

As  seen  in  Figure  15,  the  predicted  swell  arrival  time 
nearly  matches  the  time  shown  by  both  sets  of  observations. 
However,  the  predicted  energy  of  the  rise  is  slightly  higher 
than  the  observational  values.   The  peak  energy  (around  the 
time  of  2  0  GMT  August  15)  is  also  the  same  as  the  observa- 
tional peak. 

After  00  GMT  August  16,  the  prediction ' site  is  within 
the  typhoon's  wind  circulation  area.   Most  wave  heights (H  1/10) 
are  8m  and  the  highest  observation  is  10m.   These  waves 
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dominate  the  entire  wave  field  and  no  swell  can  be  dis- 
tinguished later  than  00  GMT  August  16. 
3.   Typhoon  Owen 

In  Typhoon  Owen  (see  Figure  17) ,  the  swell  arrival 
time  indicated  by  observations  is  earlier  than  the  predicted 
time  by  about  6-12  hours.   The  predicted  peak  energy  lies 
close  to  the  observational  values. 

In  fact,  the  NODC  values  are  higher  than  those  of 
the  prediction,  and  ROK  Navy  observations  are  lower. 
Except  near  the  peak  both  observation  sets  show  similar 
time  variations.   Therefore,  the  observed  values  are  con- 
sidered consistent  and  probably  accurate. 

F.   ERROR  SOURCES 

I  have  made  several  assumptions  in  this  study  in  order 
to  simplify  calculations.   The  most  serious  error  sources 
involve  assumptions  about  the  windows  in  the  Ryukyu  Islands, 
group  velocity  in  shallow  water,  and  simple  bottom  contours 
in  shallow  water. 

There  are  also  differences  between  the  predictions  (which 
include  waves  dependent  on  the  local  weather  conditions)  that 
make  it  difficult  to  evaluate  the  prediction.   Regarding  the 
assumptions  about  group  velocity  in  Chapter  III.B,  I  used 
n  =  3/4  in  shallow  water  for  calculation  of  group  velocity 
C   =  nc.   In  Typhoon  Hope,  which  has  the  longest  shallow 
water  travel  distance,  miscalculation  of  C  would  have  its 

g 
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greatest  effect.   Yet  the  observed  and  predicted  energy 
peaks  are  not  greatly  separated.   This  may  mean  that  the 
approximations  are  realistic. 

The  bottom  contour  assumptions  seem  to  be  reasonable 
through  three  tests.   The  local  weather  condition  is  the 
most  serious  factor  contributing  to  differences  between 
the  observations  and  the  predictions.   The  local  wind 
pattern  was  the  following  at  Cheju-do,  according  to  the 
ROK  Navy  data  set. 

1979   072900-073100:   SW  5-8  kts  (preceding  Hope 

swell) 

080100-080300:   SW  6-10  kts  (during  Hope 

swell) 

1979   081200-081300:   S  6-8  kts  (preceding  Irving 

swell) 

081400-081700:   S  30-60  kts  (during  Irving 

swell) 

092400-092600:   E  5-8  kts  (preceding  Owen 

swell) 

092600-093000:   E  20-30  kts  (during  Owen 

swell) 

Therefore,  before  the  swell  arrivals,  the  local  wave  heights 

were  considered  less  than  1-1.5  m.   During  Irving  and  Owen 

there  were  important  local  sea  contributions  to  observed 

height. 

Lastly,  the  original  error  sources,  i.e.,  those  involving 

input  parameters  for  the  TYWAVES  model,  the  location  of  the 

typhoon  center,  the  winds  in  the  source  region,  and  typhoon 
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size,  etc.,  are  ignored  in  this  study.   Those  errors  are 
discussed  in  Refs.  [2]  and  [17]. 
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V.  CONCLUSIONS 

From  Figures  13,  15  and  17 ,  I  can  conclude  that  most 
predicted  swell  heights  are  lower  than  those  of  the  obser- 
vations (combined  sea  and  swell) . 

The  times  of  occurrence  of  the  predicted  peak  heights 
agreed  reasonably  well  with  those  of  observations  for  the 
swell  from  each  typhoon. 

These  results  on  the  basis  of  this  limited  test  suggest 
that  TYWAVES  predict  satisfactorily  those  swells  in  the  East 
China  Sea  which  originate  in  tropical  cyclones  in  the 
western  North  Pacific. 

Computer  aided  predictions  may  improve  the  quality 
of  the  forecasts  by  reducing  the  need  for  simplifying 
approximations/  as  in  the  case  of  the  treatment  in  this 
thesis  of  the  shoaling  and  refraction  processes.   Such 
predictions  would  also  provide  a  larger  base  for  assessing 
the  accuracy  of  the  method.   Additional  verification  is 
required  to  draw  more  specific  conclusions. 
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APPENDIX  A 

Outputs  of  TYWAVES  For  Typhoon  Owen 
at  12  GMT  September  26 

The  following  tables  and  figures  show  the  computer  anal- 
ysis for  Typhoon  Owen  by  TYWAVES. 

a.   Period-directional  spectrum  at  each  of  12  sources. 

b-1.  The  significant  wave  height  (H  1/3  in  feet)  distri- 
bution in  typhoon  area.   "0"  indicates  the  land 
area  and  the  distance  between  the  grid  points  is 
40  NM. 

b-2.  The  maximum  wave  period  distribution  in  typhoon 
area.   "-1"  indicates  the  land  area. 

b-3.  The  maximum  wind  wave  directions  distribution. 

The  direction  indicated  with  16  unit  point  rose,  from 
North  (1)  to  NNE  (16)  with  CCW  direction. 

c.    NEPRF  Typhoon  Wave  Program  Analysis  for  Typhoon  Owen 
during  the  typhoon  period  from  the  time  of  the  first 
typhoon  warning  to  72  hours  later,  based  on  post- 
analysis  data. 
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APPENDIX  B 

Sea  Spectra  at  Source  Region  for  Three  Typhoons  and  Significant 

Wave  Distribution 


The  following  tables  and  figures  show  the  period-direction 
spectrum  at  selected  source  grid  point  and  NEPRF  typhoon  wave 
program  analysis  for  wave  height  (H  1/3  in  feet)  distribution 
around  the  typhoon  center,  respectively. 

a-1.   NEPRF  Typhoon  Wave  Model  (period-direction  spectrum 

at  each  selected  grid  point)  for  Typhoon  Hope. 
a-2.   NEPRF  Typhoon  Wave  Program  Analysis  for  Typhoon  Hope. 
b-1.   NEPRF  Typhoon  Wave  Model  (period-direction  spectrum 

at  each  selected  source  grid  point)  for  Typhoon 

Irving. 
b-2.   NEPRF  Typhoon  Wave  Program  Analysis  for  Typhoon 

Irving. 
c-1.   NEPRF  Typhoon  Wave  Model  (period-direction  spectrum 

at  each  selected  source  grid  point)  for  Typhoon 

Owen. 
c-2.   NEPRF  Typhoon  Wave  Program  Analysis  for  Typhoon 

Owen. 
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APPENDIX  C 
Propagation  Energy  Arrival  Time 

The  following  tables  show  the  components  from  each 
source  grid  point,  travel  time  and  arrival  time  (GMT). 
R  and  R  are  determined  from  Map  1  (U.S.  Navigation  Chart 
No.  94027,  Scale  1:  927.700  at  lat  32°15')  from  the  critical 
depth  of  each  period. 
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Jul 
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30 

PT4 
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Ro/R 
t 
tar 
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10 
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162033 

161042 

160147 

152321 

152213 

00 

Aug 

16 

PT4 
29.6-126.0 

Ro/R 

t 

218/- 
20.57 

218/- 
14.34 

-/218 

7.41 

-/218 
6.49 

-/218 
6.09 

tar 

162034 

161421 

160725 

160629 

160605 

PT8 
27.6-126. 

Ro/R 

t 

337/- 
31.79 

337/- 
22.17 

20/317 

11.85 

-/317 
10.03 

-/317 
9.41 

tar 

170747 

162210 

161151 

161002 

160925 
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PT12 
29.6-127.5 

Ro/R 

t 

tar 

221/- 
20.85 
162051 

221/- 
14.54 
161432 

-/221 

7.52 

160731 

-/221 

6.58 

160635 

-/221 

6.17 

160610 

06 

Sep 

25 

PT4 
21.3-132.4 

Cgo/Cgs 
Ro/R 
t 
tar 

;  10.6 
784 
73.96 
280758 

15.2 
784 
51.58 
270935 

18.0/29. 
654/730 
39.39 
262123 

4  24.3/33. 
621/163 
30.44 
261225 

6  28.8/35 
611/173 
26.05 
260803 

PT12 
21.3-133.9 

Ro/R 

t 
tar 

822 

77.55 

281133 

822 

54.08 

271205 

702/120 

41.62 

262337 

662/160 

32.00 

261400 

658/164 

27.43 

260926 

18 

Sep 

25 

PT4 
22.6-131.7 

Ro/R 

t 
tar 

696 

65.66 

281140 

696 

45.73 

271547 

563/133 

74.63 

270438 

513/183 

26.56 

262034 

508/188 

22.89 

261653 

PT12 
22.6-133.1 

Ro/R 

t 

tar 

731 

68.96 

281458 

731 

48.09 

271805 

610/121 

36.74 

270644 

577/154 

28.32 

262219 

574/157 

24.32 

261819 

06 

Sep 

26 

PT4 
23.5-131.4 

Ro/R 

t 

tar 

639 

60.28 

281817 

639 

42.04 

280002 

509/130 

31.64 

271338 

469/170 

24.36 

270622 

419/220 

20.07 

270204 

PT12 
23.5-132.8 

Ro/R 

t 

tar 

675 

63.68 

282141 

675 

44.41 

280225 

555/120 
33.76 

271546 

522/153 

26.03 

270802 

515/160 
22.35 

270421 

18 

Sep 

26 

PT4 
24.4-131.6 

Ro/R 

t 

tar 

542 

51.13 

282108 

542 

35.66 

280539 

417/125 

26.55 

272034 

382/160 

20.48 

271429 

371/171 

17.66 

271140 

PT12 
24.4-135.3 

Ro/R 

t 

tar 

667 

62.92 

290855 

667 

43.88 

281353 

607/60 

34.5 

280430 

562/. 05 

26.26 

272016 

562/105 

22.44 

271626 

06 

Sep 

27 

PT4 
25.5-131.9 

Ro/R 

t 

tar 

544 

51.32 

290919 

544 

35.79 

281747 

484/60 

27.92 

280955 

410/134 

20.86 

280252 

407/137 

17.96 

272358 

PT12 
25.5-133.4 

Ro/R 

t 

tar 

591 

55.75 

291345 

591 

38.88 

282053 

526/65 

30.34 

281220 

446/145 

22.67 

280440 

443/148 

19.52 

280131 

18 

Sep 

27 

PT4 
26.5-132.0 

Ro/R 

t 

tar 

496 

46.79 

291648 

496 

32.63 

290238 

430/66 

25.24 

281914 

348/148 

18.73 

281244 

345/151 

16.20 

281012 
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PT6 
28.5-132.1 

Ro/R 

t 
tar 

406 

38.30 

290818 

406 

26.71 

282043 

306/100 

19.76 

281346 

251/155 

14.94 

280857 

245/160 

13.01 

280706 

PT12 
26.5-133.5 

Ro/R 

t 
tar 

547 

51.60 

292136 

547 

35.99 

290559 

487/60 

28.08 

282205 

417/130 

21.03 

281502 

407/140 

18.04 

281202 

06 

Sep 

28 

PT4 
27.3-132.1 

Ro/R 

t 

tar 

460 

43.40 

300124 

460 

30.26 

291216 

400.60 
23.43 

290526 

318/142 

17.31 

282319 

315/145 

14.99 

282059 

PT12 
27.3-133.5 

Ro/R 

t 
tar 

511 

48.21 

300612 

511 
33.62 

291537 

414/97 

25.44 

290726 

386/125 

19.61 

290136 

381/130 

16.86 

282252 

18 

Sep 

28 

PT4 
27.8-132.1 

Ro/R 

t 

tar 

436 

41.13 

301108 

436 

28.68 

292241 

371/65 

22.05 

291603 

291/145 

16.29 

291017 

289/147 

14.14 

290608 

06 

Sep 

29 

PT4 
28.5-132.4 

Ro/R 

t 
tar 

416 
39.25 

302115 

416 

27.37 

300922 

361/55 

21.18 

300311 

291/125 

15.70 

292142 

288/128 

13.58 
291935 
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APPENDIX  D 
Shoaling  of  the  Spectral  Energy  Components 

For  refraction  factor  computation  at  00  GMT  August  13 

from  Point  4  of  Typhoon  Irving,  the  critical  water  depth  is 

32S  feet  and  the  angle  between  the  swell  crest  and  the 

bottom  contour  of  325  ft  depth  is  40°.   Thus,  with  h/gT2 

and  the  Figure  2-19  [7],  the  refraction  factors  are  derived 

like  the  following: 

T  (sec)  =7    10    13    16    19 

K  (-— )  =1     1     1   0.98  0.963 
r  h 
o 

In  the  same  way,  the  other  refraction  factors  are  shown  in 
shoaling  computation,  in  the  following  tables: 

1.  Shoaling  Computation 

a.  Typhoon  Hope 

b.  Typhoon  Irving 

c.  Typhoon  Owen 

2 .  Shoaling   Energy  and  Arrival  Times 

a.  Typhoon  Hope 

b.  Typhoon  Irving 

c .  Typhoon  Owen 
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1.  Shoaling  Conputation 


a.  Typhoon  Hope 

12  Jul  30    PT4 


T 
k 


PT12 


a  =35 
o 


00  Jul  31    PT4 


SSE. 


a  =  45     SSE 
o  s 


PT12 


(ovrr)   17.4-133.9  (k  k  r  1 

h  =  413     SSE      7 
o 

a  =40°     SSE 
o  s 


17.4-135.3  (k  k  r  1 
s  r 

h  =  413     SSE     8 


18.6-131.5   (k  k  )2  1 
s  r 

h  =  413     SSE     4 
o 


10   13 
1   0.92 


E  (2  •  kilo- 
16      19   H  1/10  (in 

0.89    0.90   ^A^ 

meters J 


1  1      0.989   0.972 

1  0.8464  0.7748  0.7653  E    H  1/10 


19  11 
19  9.3 


9 
7.0 


14     63 
10.7   53 


3.17 
2.91 


11  1  0.98 

1  0.8464  0.7921  0.7779 

28  31  36  42     147   4.85 

28  26.2  28.5  32.7   123.4  4.44 

1  1  0.984  0.962 

1  0.8464  0.7670  0.7496 

15  20  13  15     69    3.32 

15  16.9  10.0  11.2   57.1  3.02 


1  1 


0.982 


18.6-132.9   (k  k  r  1 
s  r 


h  =  431     SSE. 


a  =  35     SSE 
o 


1  0.8467  0.7921  0.7811 
37  43      19     19     123   4.44 


37  36.4    15 


14.8   106.2  4.12 


12  Jul  31    PT4 


K 


19.6-128.3   (kskrr  1 


n  =  413 


a  =  45 
o 


1  1  0.984  0.962 

1  0.8464  0.767  0.75 

23  30  46  20 

23  25.4  35.3  15 


121   4.40 
98.7  3.97 


00  Aug  01    PT4 


20.6-123.2   (k  k  )   1 
s  r 


h  =  236     SSW 


a  =  60     SSW 
o  s 


1  0.935  0.875  0.825 

1  0.7399  0.6065  0.5513 

23  30  19  40     118   4.35 

23  22.2  11.5  22.1   83.8  3.66 


95 


PT12 

k 

r 

1 

20.6-126. 

8 

(kk)2 
s  r 

1 

h  =  413 

S 
o 

0 

a    =  35 
o 

S 

s 

0 

12  Aug  01    PT4 

k 

r 

1 

21.5-122. 

2 

(kk)2 
s  r 

1 

h  =  236 

ssw 

o 

0 

a  =  30° 
o 

SSW 
s 

0 

PT8 

k 
r 

1 

19.5-122. 

2 

(kk)2 
s  r 

1 

h  =  295 

SSW 
o 

1 

rrO 

a  =  55 
o 

SSW 
s 

1 

PT12 

k 

r 

1 

21.5-123. 

7 

(k  k  )2 
s  r 

1 

h  =  295 

SSW 
o 

0 

a  =  50° 
o 

SSW 
s 

0 

b.  Typhoon  Irving 

00  Aug  13    PT4 

(GMT)    20.0-128.8 


11     1  0.98 

I  0.8464  0.7921  0.7779 

20  71     98  52     243   6.24 

20  60.1   77.6  40.5   198.2  5.63 

1  0.99   0.98  0.97 

1  0.8296  0.7607  0.7621 

II  33     3  10     58    3.05 
11  27.4   2.3  7.6    28.3  2.78 

1  0.97   0.93  0.88 

1  0.7963  0.6851  0.6272 

22  51     8  13     98    3.96 

22  40.6   5.5  8.2    77.3  3.52 

1  0.977  0.95  0.915 

1  0.8079  0.7149  0.6781 

6  17     49  31     104   4.08 

6  13.7   35.0  21     75.7  3.48 


deep       SSE  3    39  46     21  27     137   4.68m 

shoal  k  111  0.98  0.963 

r 

h=  325  ft      (k  k  )2  1         1       0.8464  0.7607  0.7512 


s  r 

=  40°  SS 

o 


a    =  40°         SSE  3       39       38.9       16.0         20.3       117.2  4.33 
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PT12 

20.0-130.2 

deep       SSE     5    25  38     20      11     100   4.00 

shoal      k       1111      1 

h  =  325  (kk)2     1  1     0.8464  0.7921     0.81 

s  r 

a    =  10  SSE  5         25     35.0       15.8         8.9         89.7     3.79 

o 


12  Aug  13         PT4 

22.0-128.2 


deep  SSE  4         30     31  25  13  105       4.10 

shoal  k  1111  1 

r 

h  =  325     (kk)2  1     1  0.8464  0.7921  0.81 
s  r 

a    =  10     SSE     4    30  26.2   19.8    10.5   90.5  3.81 
o 

PT12 

22.0-129.6 

deep       SSE     3    24  20     15      4      67    3.27 

shoal      k       11111 
r 

h  =  325  (kk)2     1  1       1  0.8464     0.79210.81 

s  r 

a    =  10  SSE  3         24     16.9       11.9         3.2         59         3.07 

o 


00  Aug  14         PT4 

23.5-125.0 


deep  &  4         17     22  18  17  79         3.56 

shoal  k  111  0.97         0.95 

r 

h  =  325  (kk)2     1  1     0.8464  0.7453     0.7310 

s  r 

a    =  4.5  S  4         17     is. 6       13.4         12.4       65.4     3.23 

o 
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PT12 

23.5-128.6 

deep       SSE 

4 

shoal      k 

1 

30  30     28  20     114   4.27 

1  1     0.98  0.963 

h  =  325     (kk)2  1     1  0.8464  0.7607  0.7512 
s  r 


a  =  40 
o 

SSE 

4 

30 

25.4 

21.3 

15.0 

95.7 

3.91 

Aug  14    PT4 

24.6-126. 

.7 

deep 

S 

1 

19 

31 

13 

8 

74 

3.44 

shoal 

k 

r 

1 

1 

0.975 

0.94 

0.91 

h  =  266 

(kk)2 
s  r 

1 

1 

0.8046 

0.7 

0.6708 

a  =  50 

o 

S 

1 

19 

24.9 

9.1 

5.4 

59.4 

3.08 

PT12 

24.6-128. 

,2 

deep 

SSE 

1 

28 

33 

34 

27 

125 

4.47 

shoal 

k 

1 

1 

1 

0.98 

0.963 

h  =  325     (kk)2  1     1  0.8464  0.7607  0.7512 


s  r 


a     =  40 
o 

SSE 

1 

00  Aug  15    PT4 

25.9-126.5 

deep 

S 

9 

shoal 

k 

1 

28  27.9   25.9    20.3   103.1  4.06 


40  36     23     22     132   4.60 

1  0.99   0.976   0.948 
r 

h  =  266     (kk)2  1     1  0.8296  0.7545  0.728 
s  r 

a  =  40     Q        9    40  29.9   17.4    16.0   112.3  4.24 
O  o 
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PT12 

25.9-128. 

,0 

deep 

SSE 
o 

3 

shoal 

k 

r 

1 

h  =  325 

(kk)2 
s  r 

1 

a  =  45 
o 

SSE 
s 

3 

deep 

S 

3 

shcal 

k 

r 

1 

h  =  472 

(kk)2 
s  r 

1 

a  =  50 
o 

S 

3 

12  Aug  15    PT4 

27.5-126. 

,0 

deep 

S 

0 

shoal 

k 
r 

1 

h  =  266 

(kk)2 
s  r 

1 

a  =  45 
o 

S 

0 

PT8 

25.5-125. 

,9 

deep 

s 

1 

shoal 

k 

1 

29  20     17  40     110   4.20 

11     0.98  0.95 

1  0.8464  0.7607  0.7310 

29  16.9   12.9  29.2   91.0  3.82 


34  41     33  12     125   4.47 

1   1     0.985  0.96 

1  0.8464  0.7685  0.7465 

34  34.7   25.4  9.0    106.1  4.12 


19  13     24  15     73    3.42 

1  0.983  0.957  0.93 

1  0.8179  0.7254  0.7006 

19  10.6   17.4  10.5   57.5  3.03 


17  15     14  13     62    3.15 

1  0.988  0.967  0.95 

r 

h  =  261     (kk)2  1     1  0.8262  0.7407  0.7310 

a  =  40    s       1    17  12.4   10.4  9.5    50.3  2.84 
o        ° 
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PT12 

29.5-127. 

.5 

deep 

SSE 

1 

shoal 

k 
r 

1 

h  =  325 

(kk)2 
s  r 

1 

a  =  45 
o 

SSE 

1 

PT12 

27.5-127. 

,5 

deep 

S 

0 

shoal 

k 

1 

21  56     55  75     212   5.82 

1  1     0.97  0.95 

1  0.8464  0.7353  0.7310 

1  47.4   41  54.8   145.2  4.82 


17  38     50  43     153   4.95 

1  1     0.975  0.94 

r 

h  =  384     (k  k  )2  1     1  0.8464  0.7530  0.7157 


s  r 


a  =  51     ^ 
o 

0 

00  Aug  16    PT4 

29.6-126.0 

deep       S 

0 

shoal      k 

1 

17     32.2       39.9         30.8       119.9  4.38 


15     14  19  15  65         3.22 

111  1 


h  =  266     (kk)2  1     1  0.8464  0.7921  0.81 


s  r 


a  =  0 

o 

S 

1 

PT8 

27.6-126.0 

deep 

s 

0 

shoal 

k 
r 

1 

h  =  266 

(kk)2 
s  r 

1 

a  =  45 
o 

s 

0 

15  11.8   15.0    12.2   54    2.94 


16  20     31  21     90    3.79 

1  0.983  0.957  0.93 

1  0.8179  0.7254  0.7006 

16  16.4   22.5  14.7   69.6  3.34 


100 


PT8 

27.6-126.0 

deep       SSW 

0 

shoal      k 

1 

33     62  47  31  175       5.29 

1     0.976     0.93         0.87 
h  =  325  (kk)2     1  1     0.8063  0.6851     0.6131 


a    =  60 
o 

SSW 

0 

33 

50 

32.2 

19 

134.2 

4.63 

PT12 

29.6-127. 

.5 

deep 

SSE 

1 

24 

54 

60 

38 

179 

5.35 

shoal 

k 

r 

1 

1 

1 

1 

1 

h  =  354 

(kk)2 
s  r 

1 

1 

0.8464 

0.7921 

0.81 

a     =  15 
o 

SSE 

1 

24 

45.7 

47.5 

30.8 

149 

4.88 

PT12 

29.6-127. 

,5 

deep 

S 

1 

21 

40 

40 

19 

122 

4.42 

shoal 

k 
r 

1 

1 

1 

0.96 

0.93 

h  =  384 

(kk)2 
s  r 

1 

1 

0.8464 

0.73 

0.7001 

a  =  50 
o 

S 

1 

21 

33.9 

29.2 

13.3 

98.4 

3.97 
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c. 

Typhoon  Owen 

06 

Sep  25 

PT4 

kr    2 
(ksV 

1 

1 

1 

1 

0.98 

(GMT) 

21.3-132. 

.4 

1 

1 

0.8464 

0.7921 

0.7779 

h  =  413  ft 

(SSE)Q 

7 

26 

31 

8 

0 

74 

3.44 

"SCO 

a    =  35 
o 

(SSE) 

7 

26 

26.2 

6.3 

0 

65.5 

3.24 

PT12 

k 
r 

1 

1 

1 

1 

0.98 

21.3-133. 

,9 

(kk)2 
s  r 

1 

1 

0.8464 

0.7921 

0.7779 

h  =  413 

SSE 
o 

5 

28 

19 

14 

10 

78 

3.53 

a  =  35 
o 

SSE 
s 

5 

28 

16.1 

11.1 

7.8 

68 

3.30 

18 

Sep  25 

PT4 

k 
r 

1 

1 

1 

0.989 

0.973 

22.6-131. 

,7 

(kk)2 
s  r 

1 

1 

0.8464 

0.7335 

0.7716 

h  =  413 

SSE 
o 

3 

38 

36 

10 

0 

89 

3.77 

a  =  40 
o 

SSE 
s 

3 

38 

30.5 

7.3 

0 

78.8 

3.55 

PT12 

k 
r 

1 

1 

1 

1 

0.985 

22.6-133. 

,1 

(kk)2 
::s  r 

1 

1 

0.8464 

0.7921 

0.7859 

h  =  413 

SSE 
o 

2 

27 

34 

15 

9 

88 

3.75 

a  =  30 

o 

SSE 

s 

2 

27 

28.8 

11.9 

7.1 

76.8 

3.50 

06 

Sep  26 

PT4 

k 
r 

1 

1 

1 

0.985 

0.95 

23.5-131. 

,4 

(kk)2 
s  r 

1 

1 

0.8464 

0.7685 

0.73 

h  =  384 

SSE 
o 

2 

39 

54 

42 

1 

140 

4.73 

a  =  47 

0 

SSE 
s 

2 

39 

45.7 

32.3 

0.7 

119.7 

4.38 
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PT12  kr  111  1  0.985 

23.5-132.8   (kskr)2  1   1   0.8464  0.7921  0.7859 

h  =  384     SSE     2   33  49     26      19      132   4.60 
o 

a  =  33     SSEo     2   33  41.5    20.6    14.9    112   4.23 
o  s 

18  Sep  26  PT4  k  111  1  0.983 

24.4-131.6   (k  k  )2  1   1   0.8464  0.7921  0.0827 

h  =  384     SSE     0   32  54     48     53     188   5.48 
o 

a  =  35     SSEo     0   32  45.7    38      41.5    157.2  5.01 
o  s 

PT12  k  111  1  1 

r 

24.4-135.3      (k  k  )2     1       1       0.8464     0.7921     0.81 
s  r 

h  =  431  SE  0       9       15  17  15 

o 

a    <10  SE  0       9       12.7         13.5         12.2         47.4     2.75 

o  s 

06  Sep  27  PT4  k  111  1  1 

25.5-131.9      (k  k  )2     1       1       0.8464     0.7921     0.81 
s  r 

h  =  354  SSE  0       17     54  31  56  150       4.90 

o 

a    <10  SSE  0       17     45.7         24.6         45.4         132.7  4.61 

o  s 

PT12  k  111  1  1 

r 

25.5-133.4      (k  k  )2     1       1       0.8464     0.7921     0.81 
s  r 

h  =  413  SE  1       9       11  7  8  37         2.43 

o 

a  <10  SE  19       9.3  5.5  6.5  31.3     2.23 

o  s 
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SSE. 


h  =  431 
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SE 


SE 


PT6 
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11  1  1 

I  0.8464  0.7921  0.81 

30  58  49  39      180   5.37 

30  49.1  38.8  31.6    150.5  4.91 

II  1  1 

1       0.8464  0.7921  0.81 

18     39  32  30  123       4.44 

18     33  25.3  24.3         101.6  4.03 
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11  1  1 
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75         3.46 
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SE. 
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a  <0 
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SE 


SE 
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11  1  1 

I  0.8464  0.7921  0.81 

26  57  81  35      201   5.67 

26  48.2  64.2  28.4    166.8  5.17 

II  1  1 

1       0.8464  0.7921  0.81 

26     44  52  36  161       5.08 

26     37.2  41.2  29.2         134.6  4.64 
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Shoaling   Energy    and  Arrival   Time 
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T 

7 
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+--.-W 
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11  1  1 

I  0.8464  0.7921  0.81 
15     25  24  14 
15     21.2  19  11.3 

II  1  1 

I  0.8464  0.7921  0.81 
10     15  16  5 

10     12.7  12.7  4.0 

II  1  1 

1       0.8464  0.7921  0.81 

10     18  12  11 

10     15.2  9.5  8.9 


10 


13 


16 


81         3.60 
67.5     3.29 


49         2.80 
40.4     2.54 


55         2.96 
45.6     2.70 
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t 
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27.5-126.0 
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13 

24 

15 

73 
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161832 

160635 
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17 

15 

14 

13 

62 
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17 
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PT12 

t 

32.54 
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53 

43 

153 

Ss 

0 

17 
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t 
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6.49 

6.09 

tar 

162034 

161421 

160725 

160629 

160605 

29.6-126.0 
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31.79 
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SSW 
o 

0 

33 

62 

47 

31 

175 

SSW 

0 

33 

50 

32.2 

19 

134.2 

109 


PT12 

t 

20.85 
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